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14.  ABSTRACT 

The  use  of  piezoelectric  wafer  active  sensor  (PWAS)  for  structural  health  monitoring  (SHM)  has  developed  into  a  promising 
tool  for  real-time  structural  damage  assessment.  The  1-year  research  work  described  in  this  report  addresses  the  following 
issues:  (i)  PWAS  durability  and  survivability;  (ii)  optimum  PWAS  design;  (iii)  steering  beam  large-area  coverage  (embedded 
ultrasonic  structural  radar  -  EUSR);  and  (iv)  probability  of  detection  (POD)  with  PWAS  probes  (an  initial  investigation). 

The  durability  and  survivability  of  the  PWAS  transducers  has  been  confirmed  under  various  environmental  exposures  for  the 
duration  covered  in  this  one-year  effort.  These  results  give  confidence  in  this  new  technology  for  installation  on  realistic 
structures.  The  report  presents  analytical  methods  and  experimental  results  that  can  be  successfully  used  in  the  design  of 
active  SHM  systems.  The  EUSR  steering  beam  PWAS  phased-arrays  method  has  been  refined  to  increase  it  power  of 
detection  for  various  structural  defects.  It  has  been  also  proved  that  EUSR  performance  is  not  affected  by  typical  structural 
curvature.  The  initial  PWAS  POD  investigation  has  indicated  that  crack  as  small  as  1.57  mm  can  be  detected  and  imaged 
relatively  easily.  Overall,  the  results  confirm  that  PWAS  transducers  are  a  technology  with  good  potential  for  implementation 
in  active  SHM  systems. 
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2.  SUMMARY 

This  report  consists  of  four  unique  tasks  that  serve  to  evaluate  the  piezoelectric  wafer  active  sensors’ 
(PWAS)  performance,  durability,  design,  and  reliability.  Task  1:  PWAS  Durability  and  Survivability 
showed  how  temperature,  weather,  liquids,  and  strain  affect  the  performance  of  the  PWAS.  Minimal 
changes  were  caused  by  temperature,  weather,  and  liquids,  but  it  was  discovered  that  large  changes  in 
impedance  can  be  caused  by  strain.  Task  2:  Optimum  PWAS  Design  uses  theoretical  model  and  real  tests 
to  design  a  PWAS  layout  that  will  detect  cracks  and  corrosion.  MathCAD  and  MatLAB  were  used  to 
create  models  of  Lambs  in  a  structure  and  experiments  were  done  to  verify  the  results.  Task  3:  Steering- 
Beam  Large  Area  Coverage  examined  what  affect  multiple  cracks,  crack  orientation,  and  curvature  has  on 
PWAS.  New  techniques  like  the  Hilbert  transformation  and  the  autoswitch  device  helped  to  improve  the 
scanning  and  detection  methods.  The  last  task.  Task  4:  Probability  of  Detection  (POD),  uses  pin  holes  to 
test  the  smallest  amount  of  damage  the  PWAS  can  detect  which  was  about  1.57  mm. 

3.  INTRODUCTION 

Structural  health  monitoring  (SHM)  is  a  method  of  determining  the  health  of  a  structure  from  the  readings 
of  an  array  of  permanently-attached  sensors  that  are  embedded  into  the  structure  and  monitored  over  time. 
SHM  can  be  performed  in  basically  two  ways,  passive  and  active.  Passive  SHM  consists  of  monitoring  a 
number  of  parameters  (loading  stress,  environment  action,  performance  indicators,  acoustic  emission 
from  cracks,  etc.)  and  inferring  the  state  of  structural  health  from  a  structural  model.  In  contrast,  active 
SHM  performs  proactive  interrogation  of  the  structure,  detects  damage,  and  determines  the  state  of 
structural  health  from  the  evaluation  of  damage  extend  and  intensity.  Both  approaches  aim  at  performing 
a  diagnosis  of  the  structural  safety  and  health,  to  be  followed  by  a  prognosis  of  the  remaining  life.  Passive 
SHM  uses  passive  sensors  which  only  “listen”  but  do  not  interact  with  the  structure.  Therefore,  they  do 
not  provide  direct  measurement  of  the  damage  presence  and  intensity.  Active  SHM  uses  active  sensors 
that  interact  with  the  structure  and  thus  determine  the  presence  or  absence  of  damage.  The  methods  used 
for  active  SHM  resemble  the  methods  of  nondestructive  evaluation  (NDE),  e.g.,  ultrasonics,  eddy 
currents,  etc.,  only  that  they  are  used  with  embedded  sensors.  Hence,  the  active  SHM  could  be  seen  as 
“embedded  NDE”.  One  widely  used  active  SHM  method  employs  piezoelectric  wafer  active  sensors 
(PWAS),  which  send  and  receive  Lamb  waves  and  determine  the  presence  of  cracks,  delaminations, 
disbonds,  and  corrosions.  Due  to  its  similarities  to  NDE  ultrasonics,  this  approach  is  also  known  as 
“embedded  ultrasonics”. 

Using  Lamb  waves  (guided  waves)  is  much  more  efficient  for  inspecting  thin- wall  structures  than  using 
conventional  pressure  waves.  Lamb  waves  remain  guided  between  two  parallel  free  surfaces,  such  as  the 
upper  and  lower  surfaces  of  a  plate  or  shell  and  can  propagate  along  the  mid-surface  of  thin-wall  plates 
and  shallow  shells.  They  can  travel  at  relatively  large  distances  with  very  little  amplitude  loss  and  offer 
the  advantage  of  large-area  coverage  with  a  minimum  of  installed  sensors.  In  Lamb-waves  applications, 
PWAS  couple  their  in-plane  motion,  excited  through  the  piezoelectric  effect,  with  the  Lamb-waves 
particle  motion  on  the  material  surface.  Thus,  PWAS  can  be  both  exciters  and  detectors  of  elastic  Lamb 
waves  traveling  in  the  material.  Damage  is  detected  if  the  exciter  and  detector  signals  are  different. 

The  electromechanical  (E/M)  impedance  method  is  an  embedded  ultrasonic  method  that  is  emerging  as  an 
effective  and  powerful  technique  for  SHM.  The  impedance  method  provides  several  advantages  over 
traditional  NDE  approaches.  Through  PWAS  permanently  attached  to  the  structure,  the  E/M  impedance 
method  is  able  to  measure  directly  the  high-frequency  local  impedance  spectrum  of  the  structure.  Because 
the  high-frequency  local  impedance  spectrum  is  much  more  sensitive  to  incipient  damage  than  the  low- 
frequency  global  impedance,  the  E/M  impedance  method  is  better  suited  for  applications  in  structural 
health  monitoring  then  other  more  conventional  methods.  The  PWAS  health  monitoring  technique  has 
been  successfully  applied  to  various  structures  ranging  from  aerospace  to  civil  structures. 
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4.  TASK  1:  PWAS  DURABILITY  AND  SURVIVABILITY 


The  purpose  of  this  task  is  to  explore  the  durability  and  survivability  issues  associated  with  various 
environmental  conditions  on  PWAS  for  impedance-based  structural  health  monitoring.  Our  approach 
consists  of  several  stages.  First,  we  put  the  PWAS  into  oven  to  determine  how  it  behaves  in  the 
temperature  range  associated  with  operational  usage.  Free  PWAS  and  bonded  PWAS  were  placed  in  an 
oven  and  exposed  to  a  certain  temperature  and  a  heating  cycle.  Second,  we  put  the  free  PWAS  and 
bonded  PWAS  with  different  adhesives  and  coatings  in  the  outdoor  environment  to  determine  how  it 
behaves  in  different  weather  conditions  (sun,  rain,  etc.).  Third,  we  determine  how  PWAS  behaves  when 
exposed  to  water  and  various  maintenance  fluids  (hydraulic  fluids  and  lubrication  oils).  Forth,  determine 
how  the  piezoelectric  material  behaves  when  exposed  to  tension  stresses  associated  with  Air  Force  usage. 

Task  Overview 

o  Environmental  Conditions  Tests 

■  Subtask  1.1-  Oven  Test 

•  Heat  cycling 

■  Subtask  1.2-  Outdoor  Test 

•  Real  World  Climate  and  Environmental  Changes 

■  Subtask  1.3  -  Submersion  Test 

•  Distilled  water 

•  Saline  water 

•  Hydraulic  fluid  MIL-PRF-  83282  Synthetic  hydrocarbon 

•  Hydraulic  fluid  MIL-PRF-  87257  Synthetic  hydrocarbon 

•  Hydraulic  fluid  MIL-PRF-  5606  Mineral 

•  Aircraft  cube  oil  MIL-PRF-7808L  Grade  3  Turbine  engine  synthetic 

•  Kerosene  Hydraulic  Fluid 
o  Physical  Properties  Tests 

■  Subtask  1.4-  Tensile  Test 


4.1.  Subtask  1.1-  Oven  Test 

4.1.1.  Methods,  Assumptions,  and  Procedures 

Oven  test  is  a  part  of  tests  of  durability  and  Survivability  of  Piezoelectric  Wafer  Active  Sensor  (PWAS). 
The  objective  of  this  test  is  to  determine  how  the  piezoelectric  material  behaves  in  the  temperature  range 
associated  with  SHM  usage.  The  test  consist  of  exposing  PWAS  in  an  oven  to  the  given  temperature 
range  for  the  specified  period  of  time.  The  functionality  of  the  PWAS  has  been  tested  before  and  after 
temperature  exposure  via  impedance  spectra  measured  using  the  Hewlett  Packard  4194 
Electromechanical  Impedance  Analyzer. 

The  impedance  of  free  PWAS  has  been  measured  on  the  following  frequency  ranges:  R1-R2:  100-900 
kHz  (in  plane  vibrations  narrow),  Rl-R#:  100-3100  kHz  (in  plane  vibrations  wide),  and  Tl:  10.5- 
12.5MHz  (thickness  vibrations).  Additionally,  PWAS  bonded  to  metallic  plate  will  be  tested  on  the  range 
P1-P6:  1-40  kHz  (plate  flexural  modes). 


Equipment/Apparatus  and  Specimens 

The  major  equipment  that  was  used  was  the  Blue  M  Electric  Oven  and  Hewlett  Packard  4194  Electro- 
Mechanical  Impedance  Analyzer  Figure  1. 
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Figure  1 


(a)  Blue  M  Electric  Oven  (Max  Temperature:  704°C),  (b)  HP  4194  Impedance  Analyzer 


Two  types  of  specimens  have  been  used  (Figure  2).  The  bond  adhesive  is  M  Bond  200. The  connection 
wire  is  AWG  38.  The  diameter  of  wire  is  0.1mm. 


Figure  2  Specimens,  (a)  Schematic  free  PWAS  and  PWAS  attached  to  metallic  plate  specimen,  (b) 
Free  PWAS  and  PWAS  attached  to  metallic  plate  specimen 


Heating  Procedure 

The  heating  tests  consist  of  two  parts:  the  first  part  is  placing  PWAS  in  an  oven  and  exposing  them  to 
175F  for  an  hour  and  cool  down  to  room  temperature  five  times.  The  other  part  is  placing  PWAS  in  an 
oven  and  exposing  them  to  an  automatic  cycling  temperature  for  a  given  amount  of  time. 
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One  Cycle  consists  of  two  segments.  Segment  1:  Rise  from  lOOF  to  175F  for  an  hour  and  last  for  5 
minutes  in  175F.  Segment  2:  Fall  from  175F  to  lOOF  for  an  hour  and  last  for  5  minutes  in  lOOF  (Figure 
3). 


Temperature  Cycle  Graph 


Figure  3  Temperature  cycles  graph 


4.1.2.  Results  and  Discussions 

4J.2,L  FreePWAS 

The  graphs  below  show  the  behavior  of  the  free  PWAS  #1.  After  a  1200-cycle  oven  exposure,  no 
significant  change  was  observed  in  the  PWAS  spectrum  (Figure  4). 


(a) 

Figure  4 
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Comparison  of  PWAS  #01  with  polyurethane  coating  spectra.  After  a  1200-cycle  oven 
exposure,  no  significant  change  is  recorded:  (a)  narrow-band  spectrum  showing  only  the  first 
two  resonance  frequencies:  (b)  broadband  spectrum  showing  several  resonance  frequencies 


4.L2.2.  Bonded  PWAS 

The  graphs  below  show  the  behavior  of  the  bonded  PWAS  #2.  After  a  1200-cycle  oven  exposure,  no 
significant  change  was  observed  in  the  PWAS  spectrum  (Figure  5). 


Figures  Comparison  of  PWAS  #02  with  polyurethane  coating  spectra.  After  a  1200-cycle  oven 
exposure,  no  significant  change  is  recorded:  (a)  narrow-band  spectrum  showing  only  the  first 
four  resonance  frequencies;  (b)  broadband  spectrum  showing  several  resonance 
frequencies 


4.2.  Subtask  1.2  -  Outdoors  Test 

4.2.1.  Methods,  Assumptions,  and  Procedures 

The  objective  of  outdoors  test  is  to  determine  how  the  piezoelectric  material  behaves  in  the  different 
weather  conditions  (sun,  rain  ...)  associated  with  SHM  usage.  The  test  consists  of  placing  Piezoelectric 
Wafer  Active  Sensors  (PWAS)  in  a  cage  and  exposing  them  to  the  outdoor  environment  in  the  parking 
lot.  The  intended  duration  of  the  test  is  one  year.  The  functionality  of  the  PWAS  was  tested  before  and 
after  environment  exposure  via  impedance  spectra  measured  using  the  Hewlett  Packard  4194 
Electromechanical  Impedance  Analyzer.  The  periodicity  of  measurements  will  be  one  per  week. 

The  impedance  of  free  PWAS  will  be  measured  on  the  following  frequency  ranges:  R1-R2:  100-900  kHz 
(in  plane  vibrations  narrow),  Rl-R#:  100-3 100  kHz  (in  plane  vibrations  wide),  and  Tl:  10.5-12.5  MHz 
(thickness  vibrations).  Additionally,  PWAS  bonded  to  metallic  plate  will  be  tested  on  the  range  P1-P6:  1- 
40  kHz  (plate  flexural  modes). 

4.2.2.  Specimens  and  Equipments 

Sixteen  specimens  were  fixed  with  a  stand  which  was  put  into  a  cage  (Figure  6).  Types  and  numbers  are 
listed  in  Table  1  and  Table  2.  The  cage  is  made  by  hand  and  the  Hewlett  Packard  4194  Electro- 
Mechanical  Impedance  Analyzer  was  used  to  record  the  data. 
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Figure  6  Fixtures  for  testing  the  PWAS  specimens:  (a)  specimens  in  a  cage  placed  outdoors;  (b) 
specimens  stand  that  fits  inside  the  cage 


Table  1  Free  PWAS 


AWG  38 

M-LINE  three  conductor 
Cables 

Protection 

coating 

No  coating 

1  (No.40) 

1  (No.42) 

M-Coat  A- 
Polyurethane 

1  (No.41) 

1  (No.43) 

M-Coat  C-Silicon 

1  (No.46) 

1  (No.44) 

M-Coat  D-Acrylic 

1  (No.48) 

1  (No.45) 

Table  2  Circular  plate  specimen  with  PWAS 


Adhesive 

M-Bond  200- 
Cyanoacrylate  w/Catalyst 

M-Bond  AE  10-2-part,  100% 
Solids  Epoxy  System 

O) 

_c 

'-4— » 

CO 

o 

o 

c 

No  coating 

1  (No.22) 

1  (No.33) 

M-Coat  A- 
Polyurethane 

1  (No.23) 

1  (No.34) 

o 

o 

CD 
-1— » 

M-Coat  C- 
Silicon 

1  (No.27) 

1  (No.35) 

o 

L_ 

D_ 

M-Coat  D- 
Acrylic 

1  (No.28) 

1  (No.36) 

Two  kinds  of  connection  wire  were  used  -  one  was  AWG  38  and  the  other  was  Micro-Measurements 
Strain  Gage  Accessories:  M-LINE  Three  Conductor  Cables.(Figure  7)  Stranded  silver-plated  copper  wire, 
3-conductor  twisted  cable,  etched  Teflon  insulation  Small,  flexible  cable.  For  use  from  -452°  to  +500°F  (- 
269°  to  +260°C)  with  color-coded  red/white/black.  The  insulation  was  treated  for  bonding. 
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Figure  7  M-LINE  Three  Conductor  Cables 
4.2.3.  Results  and  Discussions 

There  were  not  any  significant  changes  after  20-week  test  except  there  is  a  lot  of  dust  deposit  on  the 
surface  of  PWAS  (Figure  8). 


Figure  8  Samples  after  20-weeks  outdoors,  (a)  Free  PWAS,  (b)  Bonded  PWAS 


4.23,1,  Free  PWAS 

The  graphs  below  show  the  behavior  of  the  free  PWAS  #43.  After  a  20-week  outdoor  exposure,  no 
significant  change  was  observed  in  the  PWAS  spectrum  (Figure  9). 
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Comparison  of  PWAS  #43  with  polyurethane  coating  spectra.  After  a  20-week  outdoor 
exposure  no  significant  change  is  recorded:  (a)  narrow-band  spectrum  showing  only  the  first 
two  resonance  frequencies;  (b)  broadband  spectrum  showing  several  resonance  frequencies 
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4.23,2,  Bonded  PWAS 

The  graphs  below  show  the  behavior  of  the  PWAS  #23  bonded  to  a  circular  plate.  After  a  20-week 
outdoor  exposure,  no  significant  change  was  observed  in  the  PWAS  spectrum  (Figure  10). 


- Before  Outdoor 

PWAS  23  (BONDED)  — SuWoofl 

Outdoors  week 


Figure  10  Comparison  of  PWAS  #23  with  polyurethane  coating  spectra.  After  a  20-week  outdoor 
exposure  no  significant  change  is  recorded:  narrow-band  spectrum  showing  only  the  first 
four  resonance  frequencies; 


4.233.  Weather  Conditions  during  Test 


The  weather  conditions  during  the  exposure  are  recorded  from  the  NOAA  website  ( 
http://www.erh.noaa.gov/cae/climate/).  Plots  of  temperature  and  rain  over  the  testing  period  can  be 
extracted  from  these  plots  (Figure  11). 


Temperature 


MAX 

MIN 


(b) 


Figure  11  Weather  is  recorded  from  NOAA  website,  plots  from  these  data:  (a)  Temperature,  showing 
high,  low  and  average  temperature  of  days  from  initial  measurement:  (b)  Rainfall,  showing 
rainfall  of  days  from  initial  measurement 
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4.3. 


Subtask  1.3-  Submersion  Test 


4.3.1.  Methods,  Assumptions,  and  Procedures 

The  objective  of  this  test  is  to  determine  how  the  piezoelectric  material  behaves  when  exposed  to  various 
liquid  substances  associated  with  SHM  usage.  The  test  consists  of  placing  PWAS  in  containers  and 
submerging  them  in  the  substance  for  a  specified  period  of  time.  The  PWAS  are  tested  intermittently  to 
monitor  effects  from  submersion  in  the  various  fluids  (Figure  12).  Solvent  from  left  to  right:  (0)  Distilled 
water;  (1)  Saline  water;  (2)  Hydraulic  fluid  MIL-PRF-  83282  Synthetic  hydrocarbon  ;(3)  Hydraulic  fluid 
MIL-PRF-  87257  Synthetic  hydrocarbon;  (4)  Hydraulic  fluid  MIL-PRF-  5606  Mineral;  (5)Aircraft  cube 
oil  MIL-PRF-7808L  Grade  3  Turbine  engine  synthetic;  (6)Kerosene.  Impedance  spectra  acquired  over  85 
days  for  each  PWAS  are  shown  in  Figure  11,  below.  Note  that  the  impedance  spectra  of  PWAS  01, 
submerged  in  saline  water,  show  more  change  than  the  other  PWAS  spectra.  Generally,  the  spectra  for 
PWAS  01  show  the  only  significant  change  over  85  days  (Figure  13). 


Figure  12  Submersion  test  samples. 


4.3.2.  Results  and  Discussions 

Figure  13  shows  the  PWAS  impedance  spectra  collected  over  85  days  of  submergence  in  various  fluids 
associated  with  Air  Force  usage.  There  are  not  any  significant  changes  in  the  plots  for  the  synthetic 
hydraulic  fluid  and  kerosene.  On  the  saline  water  plot,  there  is  a  significant  change  on  day  85. 
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Figure  13  PWAS  impedance  spectra  coiiected  over  85  days  of  submergence  in  various  fluids 
associated  with  Air  Force  usage 
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4.4. 


Subtask  1.4-  Tensile  Tests 


4.4.1.  Methods,  Assumptions,  and  Procedures 

The  initial  step  was  to  perform  a  tensile  test  to  failure  of  one  specimen  without  sensors  in  order  to  develop 
a  stress-strain  curve  for  the  material.  The  specimen  was  placed  in  a  MTS  810  Material  Testing  System,  as 
shown  in  Figure  14  and  a  clip-on  displacement  transducer  was  mounted  to  the  specimen  (1  mm  2024  AL) 
to  measure  the  displacement. 


Figure  14  Specimen  Mounted  in  the  MTS  810  Material  Testing  System 


The  stress-strain  curve  shown  in  Figure  15  compares  the  percent  strain  to  the  stress  in  MPa. 


Stress-Strain  Curve 


Figure  15  Tensile  Specimen  Stress-Strain  Curve 
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From  the  test  data  the  material  properties  were  determined.  The  Young’s  Modulus  (£)  was  given  by  the 
slope  of  the  straight  line  region  of  the  graph  were  stress  is  proportional  to  strain,  yield  stress  (7)  was 
determined  by  the  0.2%  offset  method,  and  the  ultimate  tensile  stress  (cr^)  is  the  maximum  stress  attained 
in  the  engineering  stress-strain  diagram.  These  values  were  determined  to 

be£  «  12.5GPaJ  andcr^  «  AllMPa . 

Two  AL  specimens  were  then  tested  with  one  PWAS  attached  to  them  to  determine  how  the  PWAS  and 
the  structure  interact  under  stress.  PWAS  01  was  attached  to  the  first  AL  specimen  and  PWAS  02  was 
attached  to  the  second. 

4.4.2.  Results  and  Discussions 

4.4.2. 1.  Experimental  Results  for  PWAS  01 

PWAS  number  01  was  bonded  using  M-bond  200  cyanoacrylate  with  catalyst.  The  specimen  was 
installed  into  the  testing  machine  without  any  load  applied  and  the  baseline  impedance  was  recorded.  The 
real  part  of  the  baseline  impedance  for  the  frequency  range  of  100-200  kHz  is  shown  in  Figure  16.  The 
specimen  was  then  loaded  in  tension  until  a  strain  of  200  ps  was  obtained  and  the  new  impedance  was 
then  recorded.  This  process  was  continued  at  intervals  of  200  ps  until  5000  ps  was  achieved. 


Re(z)  vs.  Frequency 


Figure  16  Baseline  Impedance  for  PWAS  01 


Figure  17  shows  a  comparison  of  the  impedance  measurements  from  0  to  SOOOps.  Minimal  changes 
occurred  in  that  region.  Figure  18  shows  the  comparison  from  3000  to  5000ps.  Significant  changes 
occurred  above  SOOOpe.  The  specimen  was  then  unloaded  and  the  new  free  impedance  was  measured  and 
compared  to  the  initial  baseline  impedance.  This  comparison  is  shown  in  Figure  19.  It  can  be  seen  that 
after  being  loaded  to  a  strain  level  of  5000  ps  and  then  returned  to  the  unloaded  state,  the  impedance 
signature  does  not  return  to  its  baseline  impedance  signature. 
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Re(z)  vs  Frequency 
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Figure  17  Impedance  Signature  Comparisons  for  PWAS  01  for  0  ps  to  3000  pe 


Re(z)  vs  Frequency 
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Figure  18  Impedance  Signature  Comparisons  for  PWAS  01  for  3000  ps  to  5000  pc 
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Re(z)  vs  Frequency 


Figure  19  Comparison  of  Initial  Baseline  Impedance  to  Final  Unloaded  Impedance  for  PWAS  01 
4.4.2.2.  Experimental  Results  for  PWAS  02 

PWAS  number  02  was  bonded  using  M-bond  200  cyanoacrylate  with  catalyst.  Similar  to  the  specimen 
for  PWAS  01  this  specimen  was  installed  in  the  testing  machine  and  the  baseline  impedance  was  recorded 
under  zero  tension  load  condition.  The  specimen  was  then  loaded  and  impedance  readings  were  made  at 
200  |Li8  intervals  until  failure  of  the  PWAS  occurred.  The  baseline  impedance  for  PWAS  02  for  a 
frequency  range  of  100  to  200  kHz  is  shown  in  Figure  20. 
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Figure  20  Baseline  Impedance  for  PWAS  02 
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The  specimen  was  then  loaded  and  new  readings  were  taken  at  intervals  of  200|lis.  Figure  21  compares 
the  impedance  readings  for  the  baseline  up  to  3000|li8.  As  with  PWAS  01,  no  significant  changes 
occurred  up  to  the  3000|li8  level  with  exception  of  the  peak  amplitude  decreasing  with  increasing  strain. 
The  comparison  from  3000|li8  up  to  6000|li8  can  be  seen  in  Figure  22.  Significant  changes  begin  to 
happen  after  3000|li8  which  was  also  the  case  for  PWAS  01.  The  peak  amplitude  continually  drops  with 
increased  strain.  Figure  23  compares  impedance  readings  from  6000|li8  up  to  7200|li8.  The  last  reading 
that  could  be  taken  from  the  PWAS  was  at  7200|li8.  Failure  of  the  PWAS  occurred  slightly  above  that 
strain  level  and  no  further  measurements  could  be  recorded. 
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Figure  21  Impedance  Signature  Comparisons  for  PWAS  02  for  0  pe  to  3000  ps 
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Figure  22  Impedance  Signature  Comparisons  for  PWAS  02  for  3000  ps  to  6000  pe 
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Re(z)  vs  Frequency 
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Figure  23  Impedance  Signature  Comparisons  for  PWAS  02  for  6000  ps  to  7200  ps 


The  specimen  was  unloaded  after  every  lOOOjLis  and  a  new  unloaded  impedance  reading  was  recorded. 
The  new  unloaded  impedance  was  compared  to  the  initial  baseline  impedance  to  determine  if  any 
permanent  changes  had  occurred.  These  comparisons  can  be  seen  in  Figure  24  for  strain  loadings  up  to 
4000|lis  and  Figure  25  for  strains  up  to  the  failure  point. 

Figure  26  shows  a  comparison  of  the  imaginary  part  of  the  impedance  from  the  original  baseline  up  to 
4000|li8.  The  plot  continually  flattens  and  converges  to  the  smooth  curve  shown  at  4000|lis.  After  that 
point,  the  entire  curve  drops  in  value  which  can  be  seen  in  Figure  27.  Beyond  the  level  of  5000|lis  the  plot 
remains  relatively  unchanged. 
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Figure  24  Comparison  of  Initial  Baseline  to  Unloaded  Impedance  up  to  4000|a£  for  PWAS  02 
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Re(z)  vs  Frequency 
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Figure  25  Comparison  of  Unloaded  Impedance  from  4000  up  to  ZOOOps  for  PWAS  02 
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Figure  26  Imaginary  Impedance  Comparisons  for  PWAS  02  for  0  |li8  to  4000  |li8 
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Im(z)  vs  Frequency 
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Figure  27  Imaginary  Impedance  Comparisons  for  PWAS  02  for  4000  ps  to  7000  ps 


4.5.  Task  1  -  Conclusions 

The  electromechanical  (E/M)  impedance  spectra  of  free  and  bonded  PWAS  were  compared  after 
exposure  to  various  environmental  conditions,  temperature  cycles  (oven  tests),  climatic  factors  (outdoors 
tests),  operational  fluids  (immersion  tests),  and  stresses  were  used.  Various  adhesives  and  protective 
coatings  were  used  in  the  experiments. 

The  environment  tests  (Subtask  1.1  to  1.3)  were  started  about  six  months  ago  and  are  still  in  progress;  so 
far  no  significant  change  is  recorded.  This  means  the  PWAS  and  bond  between  PWAS  and  aluminum 
substrate  can  withstand  the  temperature  range  and  weather  conditions  associated  with  SHM  usage.  The 
test  will  continue  for  one  year  so  that  more  data  can  be  collected  and  analyzed.  From  the  current  results, 
PWAS  and  the  adhesive  bond  is  good  for  normal  use.  Future  environmental  condition  tests  such  as 
freeze-thaw  test  and  humidity  test  will  be  done  in  the  future. 

The  physical  properties  test  (Subtask  1 .4)  showed  that  increasing  the  strain  on  a  structure  with  the  PWAS 
attached  can  change  the  amplitude  of  the  peaks.  As  the  strain  increased,  the  peak’s  amplitude  decreased. 
Also  if  was  discovered  that  strain  alters  the  peaks  permanently  even  after  the  strain  has  stopped  being 
applied.  Strain  at  high  levels  can  also  damage  the  PWAS  and  cause  cracking.  Future  physical  properties 
tests,  like  the  shake  test,  will  be  planned  for  in  the  future. 
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5.  TASK  2:  OPTIMUM  PWAS  DESIGN 


Lamb  waves,  which  are  also  referred  to  as  guided  plate  waves,  are  a  type  of  ultrasonic  waves  that  remain 
guided  between  two  parallel  free  surfaces,  such  as  the  upper  and  lower  surfaces  of  a  plate  or  shell.  An 
alternating  voltage  applied  to  the  PWAS  terminals  produces  an  oscillatory  expansion  and  contraction  of 
the  PWAS  material.  Vice-versa,  an  oscillatory  expansion  and  contraction  of  the  PWAS  material  produces 
an  alternating  voltage  at  the  PWAS  terminals.  In  Lamb-waves  applications,  PWAS  couple  their  in-plane 
motion,  excited  through  the  piezoelectric  effect,  with  the  Lamb-waves  particle  motion  on  the  material 
surface.  Thus,  PWAS  can  be  both  exciters  and  detectors  of  elastic  Lamb  waves  traveling  in  the  material. 


Consider  a  uniform  plate  with  a  PWAS  transducer  bonded  to  its  top  surface.  When  the  PWAS  is  excited 
by  a  harmonic  electric  voltage  V  =  ,  elastic  Lamb  waves  are  excited  inside  the  structure  (Figure  28). 

The  interaction  between  the  PWAS  and  the  structure  takes  place  through  the  interfacial  shear 
stress  T  =  Tq  .  The  exact  spatial  distribution  of  Tq  (x)  depends  on  the  properties  of  the  bonding  layer 
between  the  PWAS  and  the  structure. 
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Figure  28  PWAS  interaction  with  Lamb  modes;  (a)  symmetric  Lamb  mode  So;  (b)  anti-symmetric  Lamb 
mode  Ao 


The  purpose  of  Task  2  is  the  find  the  optimum  design  for  the  PWAS.  PWAS  can  be  used  in  many  ways 
for  structural  health  monitoring  and  non-destructive  evaluation,  but  if  the  design  of  the  PWAS  is  not 
optimized,  the  effectiveness  and  performance  of  the  PWAS  can  be  compromised.  It  is  important  to  study 
the  waveforms,  surface  displacement,  material  properties  (thickness,  size,  etc),  effect  of  frequency,  effect 
the  bond,  and  even  the  sensor  shape  to  create  the  ideal  monitoring  design.  Modeling  using  MATLAB  and 
MathCAD  was  used  to  simulate  the  Lamb  waves  and  predict  their  behavior.  This  modeling  was  used  to 
find  the  optimum  PWAS  design  for  various  applications.  The  design  and  test  for  the  reciprocity  of  PWAS 
lamb  wave  transmission  and  reception  will  also  be  explored. 

Task  Overview 

■  Subtask  2.1-  Development  of  Software  Packages  for  Lamb  Wave  Simulation  in  Various  Material 
Systems 

■  Subtask  2.2  -  Lamb-Wave  Modes  Tuning  with  PWAS  Transducers 

■  Subtask  2.3  -  Modeling  and  Design  for  Reciprocity  of  PWAS  Lamb  Wave  Transmission- 
Reception 

■  Subtask  2.4  -  Modeling  and  Simulation  for  Detection  of  Corrosion  Using  Lamb  Waves 
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5.1.  Subtask  2.1  -Development  of  Software  Packages  for  Lamb  Wave  Simulation  in 

Various  Material  Systems 

5.1.1.  Methods,  Assumptions,  and  Procedures 

Software  packages  were  developed  using  MATLAB  and  MathCAD  applications.  The  software  packages 
were  intended  for  the  study  of  the  Lamb  waves  behavior  as  function  of  material  properties  and  structural 
geometry.  The  assumptions  adopted  in  this  work  were  consistent  with  the  standard  assumptions  of  the 
plain-strain  solution  to  the  wave  equation  in  a  plate  that  generate  guided  Lamb  waves.  A  full  description 
of  the  basic  assumptions  and  methods  underlying  this  approach  is  given  in  Bao  (2003). 

The  procedure  used  in  programming  these  software  packages  is  consistent  with  the  MATLAB  6.5 
programming  environment  with  graphical  user  interface  (GUI)  visual  tools. 

5. 1.1.1.  Lamb  Wave  Speed  and  Mode  Shape  Simulators 

Theoretical  modeling  of  the  equations  was  done  using  MATLAB  and  MathCAD  simulations.  Figure  29 
the  simulation  program  for  calculating  the  Lamb  wave  speed,  mode  shape,  and  stress  distribution  across 
the  thickness  and  along  a  wave  length.  Figure  30  shows  the  simulation  program  for  calculating  the  wave 
speed  and  group  velocity  for  various  material  systems  and  plate  thickness  values.  The  GUI  windows 
indicate  the  graphing  and  inputting  variables  for  Lamb  wave  simulations.  The  material  type;  Lamb  wave 
mode  type;  wall  thickness;  and  wave  speed  and  group  velocity  listed  at  selected  frequencies  can  be 
selected  in  the  windows. 


Figure  29  -  Lamb  mode  wave  speed,  mode  shape,  and  stress  distribution  simulation  program 
(MODE_SHAPE_V2) 
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Figure  30  -  (Lamb  wave  speed  and  group  velocity  simulation  program  (WAVESCOPE). 


5. 1. 1. 2.  Signal  Analysis  Software  Programs 

MATLAB  was  also  used  to  create  a  program  named  “STFT”  for  performing  signal  processing  with  the 
Short  Time  Fourier  Transform  method  (Figure  31)  and  another  software  program  named  “WAVELETS” 
for  performing  signal  processing  with  the  Continuous  Wavelet  Transform  method  (Figure  32). 
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Figure  31  -  Software  simulation  programs  for  signal  analysis  using  short  time  Fourier  transform  (STFT) 


Figure  32  -  Software  simulation  programs  for  signal  analysis  using  wavelets  (WAVELETS). 
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5. 7. 7. 3.  Math  CAD  Analysis  Software  Programs 

A  MathCAD  program  was  developed  to  predict  the  propagation  of  waves  of  various  Lamb  modes  in  a 
plate.  This  program  was  based  on  the  time-domain  Fourier  transform.  It  took  the  Fourier  transform  of  the 
excitation  signal  and  multiplied  it  with  the  Fourier  transform  of  the  wave  propagation  function  at  various 
X  locations  in  the  plate.  The  result  of  this  multiplication  was  subjected  to  the  inverse  Fourier  transform  in 
order  to  obtain  the  wave  response  at  various  x  locations  in  the  plate. 

5.1.2.  Results  and  Discussions 

The  Lamb-wave  simulation  programs  were  very  useful  for  predicting  the  behavior  of  Lamb  waves  for 
various  material  and  geometric  parameters.  Various  Lamb  wave  modes  could  be  selected  and  visualized. 
The  distribution  of  various  stresses  across  the  thickness  and  along  a  wavelength  could  be  examined.  The 
signal  processing  programs  were  found  very  useful  in  processing  the  experimental  data  and  extracting  the 
useful  information  from  noise.  The  wave  propagation  MathCAD  programs  were  found  valuable  in 
predicting  the  changes  that  appear  in  the  wave  signature  as  it  propagates  along  the  plate.  These  programs 
were  used  in  the  subsequent  subtasks  to  verify  the  optimum  sensor  design  parameters  for  various 
excitation  frequencies  and  to  understand  the  experimental  data.  All  the  above  mentioned  software 
simulation  programs  were  transmitted  to  the  AFRL  personnel  for  further  use  in  AFRL  work. 


5.2  Subtask 2.2  Lamb- Wa  ve  Modes  Tuning  with PWAS  Transducers 

In  this  section,  we  will  analyze  the  case  of  a  surface-mounted  PWAS  and  its  interaction  with  the  structural 
substrate  (Figure  33).  Because  the  PWAS  is  mounted  asymmetric  with  respect  to  the  mid-thickness  axis, 
both  symmetric  and  antisymmetric  Lamb-wave  modes  can  be  excited  and  detected.  It  will  be  shown  that 
conditions  can  be  found  for  which  the  PWAS  will  couple  preferentially  with  only  one  of  these  Lamb- 
wave  modes,  i.e.,  the  PWAS  are  tuned  with  that  particular  Lamb-wave  mode.  The  modeling  of  this 
phenomenon  is  done  as  follows. 


Figure  33  Shear  layer  interaction  between  the  PWAS  and  the  structure:  (a)  micrograph;  (b)  modeling 


5.2.1.  Methods,  Assumptions,  and  Procedures 

Consider  the  surface  mounted  PWAS  shown  in  Figure  33b.  The  PWAS  is  excited  electrically  with  a  time- 
harmonic  voltage  Ve'^"^.  As  a  result,  the  PWAS  expands  and  contracts,  and  a  time-harmonic  interfacial 
shear  stress,  develops  between  the  PWAS  to  the  structure.  This  shear  stress  acts  on  the  structure 
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over  the  interval x  g  {-a,+a) .  Thus,  the  boundary  condition  on  the  structure  upper  surface  (y  =  /z )  is 
given  by: 

=  +  (1) 

Where  H{x)  is  the  Heaviside  step  function  and  the  time-harmonic  variation  e'^^  is  implied.  From  this 

point  onwards,  the  analysis  is  performed  in  the  wave  number  domain  by  applying  the  space-domain 
Fourier  transform: 


1 

In 


(2) 


Where  —  is  the  wave  number;  the  space-domain  Fourier  transform  of  the  excitation  (1)  is 

c 


=  I  Tq{x)\^H{x  +  a)  -  H{x  -  a)^e  ^^^dx 


(3) 


Applying  the  space-domain  Fourier  transform  to  the  wave  equations  and  the  displacements,  stresses  and 
strains  we  obtain: 
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where  p  and  q  are;  p^  -  ^ 


CO 


-  ,  (j)  and  y/  are  two  potential  functions,  Cp  =  (>1  +  2//)//? 


and  cl  =  lul  p  are  the  pressure  (longitudinal)  and  shear  (transverse)  wave  speeds,  X  and  p  are  the  Lame 
constants,  p  is  the  mass  density.  The  time  dependence  is  assumed  harmonic  in  the  form  e'^^. 


Equation  (4)  accepts  the  general  solution 

sin  py  +  cos  py 

y/  =  B^  sin  qy  +  B2  cos  qy 

As  before,  the  constants  A2  and  Bi  correspond  to  symmetric  motion,  while  Ai  and  B2  correspond  to 
antisymmetric  motion.  The  excitation  can  be  also  split  into  symmetric  and  antisymmetric  components 
(Figure  34),  i.e.. 
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Figure  34  Symmetric  and  antisymmetric  loading  from  a  PWAS  of  width  2a  mounted  on  the  plate  upper 
surface:  (a)  symmetric  loading  and  (b)  anti-symmetric  loading 
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Note  that  the  negative  sign  on  the  symmetric  part  of  f  ^  ^  is  due  to  the  sign  convention  defining  the  shear 


stresses  on  the  upper  and  lower  surfaces. 


5.2.2.  Results  and  Discussions 

5.2.2. 1.  Symmetric  Solution 

The  symmetric  boundary  conditions  are: 

Because  the  motion  is  assumed  symmetric,  only  tf 
upon  substitution  into  (5),  obtain 


constants  A2  and  B\  are  retained  in  Equation  (6)  and, 


py  +  ^^1 

Uy  =  PA2  sin  py  -  sin  qy 

^yx  =  IX  [-22-^M  sin  py  +  (y  -  )B^  sin  qy'^ 

fyy  -q^  )^2  cos  py  +  li^qB^  cos 

=  -y  A2  cos  py  +  i^qB^  cos,  qy 

Substitution  of  Equations  (9)  and  (10)  into  the  boundary  conditions  (8)  yields  the  linear  system 


(9) 

(10) 

(11) 
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(12) 


(-2/^/?sin pd)A2  -q^)smqd~^B^  =~“^a 

2p 

-q^)cos pd~^A2  +{2i^qcosqd^B^  =0 


Solution  of  Equation  (12)  yields 


A, 


B, 


N, 


2p 


2p  Ds 


=  2i^q  cos  qd,  -q^)  cos  pd 

(13) 

=  (^^  -  q^  )^  cos  pd  sin  qd  +  4^^pq  sin  pd  cos 


Note  that  the  denominator  Ds  is  the  same  as  the  left  hand  side  of  the  characteristic  equation  for  symmetric 
motion.  This  indicates  that  the  poles  of  Equation  (13)  correspond  to  the  symmetric  Lamb-wave 
eigenvalues.  Substitution  of  Equation  (13)  into  Equation  (11)  yields  the  wave  number-domain 
displacement  and  strain  at  the  plate  upper  surface 


1;^=^  ^  2p  Ds 
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Where: 


Ns  =  )  cos  pd  cos  qd 


(15) 


5.2.2.2.  Antisymmetric  Solution 

The  antisymmetric  boundary  conditions  are 

i, =0 

Because  the  motion  is  assumed  antisymmetric,  only  the  constants  Ai  and  B2  are  retained  in  Equation  (6) 
and,  upon  substitution  into  (5),  obtain 


=  i^A^  sin  py  -  qB2  sin  qy 
Uy  =  pA^  cos  py  -  i^B2  cos  qy 

V  =  \}i^PA  cos  py  +  (y  -  )B^  cos  ^7] 

^yy=P  -q^)A  sin  py  +  sin  ^7] 


(17) 

(18) 


4  =  -^^A  sin  py  -  i^qB2  sin  qy  (19) 

Substitution  of  Equations  (17)  and  (18)  into  the  boundary  conditions  (16)  yields  the  linear  system 


(2/^/?sin pd^A^  +[(^^  -q^)cosqd~^B2  --^^a 
|^(^^  -^^)sin pd^A^  +{2i^qsmqd^B2  =  0 


(20) 
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Solution  of  Equation  (20)  yields 


A,=^ 


N, 


‘  2//  D 


T  N. 


2/u 


N ^  =  2i^q  sin  qd,  -^^)sin  pd 

(21) 

=  {^^  -  q^)^  sin  pd  cos  qd  +  pq  cos  pd  sin  qd 


Note  that  the  denominator  Da  is  the  same  as  the  left  hand  side  of  the  characteristic  equation  antisymmetric 
motion.  This  indicates  that  the  poles  of  Equation  (21)  correspond  to  the  antisymmetric  Lamb-wave 
eigenvalues.  Substitution  of  Equation  (21)  into  Equation  (19)  yields  the  wave  number-domain 
displacement  and  strain  at  the  plate  upper  surface 
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Where: 


Na  -  +q^)smpd  sin  qd 


(23) 


5. 2, 2, 3.  Complete  Response 


The  complete  response  to  the  PWAS  excitation  is  obtained  by  combining  the  symmetric  and 
antisymmetric  responses,  i.e., 
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We  now  return  from  the  wave  number  domain  into  the  space  domain  by  applying  the  inverse  Fourier 
transform,  i.e.. 
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Though  not  explicitly  shown,  the  functions  Ns,  Ds,  Na,  Da,  also  depend  on  co.  The  integral  in  Equation 
(26)  is  singular  at  the  roots  of  Ds  and  Da,  which  are  the  symmetric  and  anti-symmetric  eigenvalues  of  the 

Rayleigh-Lamb  equation,  i.e.,  and  Since  the  functions  Ds  and  Z)^  depend  on 

CO,  the  number  of  eigenvalues  that  exist  for  a  given  co  will  also  vary.  At  low  frequencies,  i.e.,  to^O,  only 
two  eigenvalues  exist,  and  .  At  higher  frequencies,  several  eigenvalues  will  exist.  The  evaluation  of 
the  integral  in  Equation  (26)  is  done  by  the  residue  theorem,  using  a  contour  consisting  of  a  semicircle  in 
the  upper  half  of  the  complex  plane  and  the  real  axis.  Hence: 
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The  summations  are  taken  over  all  the  symmetric  and  antisymmetric  eigenvalues  that  exist  at  the 
particular  frequency  at  which  the  analysis  is  conducted.  Similarly, 
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5. 2. 2. 4.  Closed-form  Solution  for  Ideal  Bonding 

For  thin  and  stiff  bonding  layers  ideal  bonding  between  the  PWAS  and  the  structure  may  be 
asymptotically  achieved.  For  ideal  bonding  conditions,  the  shear  stress  in  the  bonding  layer  takes  the 
simple  form  of  t^{x)  =  ar^  \^5{x- a)- S{x  +  a)^ ,  which  admits  a  closed- form  Fourier  transform,  i.e., 

T{x)y^^ I  =  ar^  \^5{x- a)- S{x  +  a)^,  f  =  ar^  [-2/ sin (29) 

Hence,  the  strain- wave  solution  becomes: 
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Similarly,  the  displacement-wave  solution  becomes: 
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(31) 


These  expressions  for  Lamb  wave  response  under  PWAS  excitation  have  been  first  published  by 
Giurgiutiu,  V.;  Bao,  J.;  Zhao,  W.  (2003).  Equations  (30)  and  (31)  contain  the  sin(^a  function.  Thus,  mode 
tuning  is  possible  through  the  maxima  and  minima  of  the  sin^^a  function.  Maxima  of  sin^^a  occur 
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when^^^j!  =  (2n-l)— .  Since ^  =  — ,  maxima  will  occur  when  the  PWAS  length  f  =2^3  equals  on  odd 

2  Z 

X 

multiple  of  the  half  wavelength  —  .  This  is  wavelength  tuning.  In  the  same  time,  minima  of  will 


occur  when  =  nji ,  i.e.,  when  the  PWAS  wavelength  is  a  multiple  of  the  wavelength.  Since  each  Lamb- 
wave  mode  has  a  different  wavespeed  and  wavelength,  such  matching  between  the  PWAS  length  and  the 
wavelength  multiples  and  submultiples  will  happen  at  different  frequencies  for  different  Lamb  modes.  An 
extensive  discussion  of  this  effect  in  terms  of  wave  number  spectra  was  given  by  Sonti  (1995).  Another 
factor  that  must  be  considered  in  Lamb  wave  tuning  under  PWAS  excitation  is  the  mode  amplitude 
tuning,  i.e.,  the  matching  of  the  PWAS  location  across  the  structural  thickness  with  the  modal 
displacement  at  that  location.  This  factor  is  contained  in  the  values  taken  by  the  function  A/D' .  For 
example,  if  the  PWAS  is  surface  mounted,  it  is  conceivable  that,  at  a  given  frequency,  some  modes  may 
have  less  surface  amplitudes,  while  other  may  have  larger  surface  amplitudes.  Thus,  two  important  factors 
for  the  design  of  PWAS-based  Lamb-wave  structural  health  monitoring  systems  has  been  identified: 
a)  The  variation  of  |  sin  \  with  frequency  for  each  Lamb  wave  mode 


b)  The  variation  of  the  surface  strain  with  frequency  for  each  Lamb  wave  mode 
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Figure  35  -  (a)  Predicted  Lamb-wave  response  a  1-mm  aluminum  plate  under  a  7-mm  PWAS  excitation: 
(a)  strain  response;  (b)  displacement  response 


A  plot  of  Equations  (30)  and  (31)  in  the  frequency  range  up  to  1000  kHz  is  presented  in  Figure  35.  (The 
So  mode  response  is  drawn  with  continuous  line,  while  the  Aq  mode  response  is  drawn  with  dotted  line.)  It 
is  apparent  that  the  strain  response  is  stronger  at  higher  frequencies  than  the  displacement  response.  This 
observation  is  consistent  with  the  fact  that  the  strain  varies  like  sin  cr,  while  the  displacement  varies 
sin  a 

like  sinccr  = - .  This  indicates  that  the  PWAS  transducers,  which  are  strain  coupled,  may  have  a  better 

a 

high  frequencies  response  than  displacement  and  velocity  transducers. 

The  Lamb-wave  mode  tuning  principles  derived  here  for  an  exciter  PWAS  apply  equally  well  for  a 
receiver  PWAS.  In  this  case,  the  receiver  PWAS  would  be  able  to  selectively  detect  certain  Lamb-wave 
modes  from  a  multimodal  Lamb  wave  reception.  An  extension  of  this  approach  to  circular  crested  Lamb 
waves  was  recently  presented  by  Raghavan  and  Cesnik  (2004)^^ 

A  MathCAD  program  for  calculating  various  axial  and  flexural  wave  parameters  was  also  created. 
(Flexural  waves  are  the  low  frequency  asymptote  of  Aq  Lamb  waves.)  The  program  calculates  the  So  and 
Ao  mode  tuning  frequency  as  function  of  plate  thickness  and  PWAS  radius.  It  also  calculates  wave  speed 
and  wave  length  as  function  of  frequency  and  plate  thickness. 


53.  Subtask  23  -Modeling  and  Design  for  Reciprocity  of  PWAS  Lamb  Wave  Transmission- 

Reception 

5.3.1.  Methods,  Assumptions,  and  Procedures 

The  objective  of  Subtask  2.3  is  to  explore  the  reciprocity  (relation  of  mutual  dependence,  action, 
interchange,  or  influence)  and  design  for  embedding  piezoelectric  wafer  active  sensors  (PWAS)  to  excite 
and  detect  tuned  Lamb  waves  for  structural  health  monitoring.  If  PWAS  can  effectively  generate  and 
detect  structural  waves,  PWAS  layouts  can  be  designed  in  arrays  to  detect  damage  by  sensing  a  change  in 
the  wave  from  one  sensor  to  the  other.  The  two  Lamb  wave  modes  that  were  used  were  symmetric  mode 
(So)  and  antisymmetric  mode  (Aq).  Before  going  into  the  procedure  section,  a  brief  review  of  Lamb- 
waves  modes  with  PWAS  theory  will  be  presented. 


5. 3. 1. 1.  Prediction  of  Reciprocity 

Theoretical  modeling  was  performed  and  then  verifled  with  experimental  tests.  A  1-D  Lamb  wave 
propagation  was  created  using  modeling  and  experimental  testing.  The  specimen  was  a  2024  T3 
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aluminum  plate  with  the  thiekness  of  1.6  mm,  PWAS  separation  of  250  mm,  and  the  exeitation  of  3-eount 
sine  tone  burst  at  a  frequeney  of  78  kHz.  Figure  36  shows  a  good  eomparison  between  the  modeling  and 
experiment. 


Figure  36  -  Plots  to  the  modeling  (a)  and  experiment  (b). 


5.3. 1.2.  Test  of  Reciprocity 


A  1.6mm  thiek,  2024  aluminum  alloy  plate  instrumented  with  an  array  of  7  mm  x  7  mm  PWAS  sensors 
(Figure  37)  was  used.  An  HP33120A  arbitrary  signal  generator  was  used  to  generate  exeitation  signals.  A 
Tektronix  TDS210  two-ehannel  digital  oseilloseope  was  used  to  eolleet  the  response  signals  from  the 
PWAS  wafer  aetive  sensors.  The  instrumentation  setup  is  presented  in  Figure  37.  The  PWAS  loeations  on 
the  plate  are  given  in  Table  3. 


Figure  37  The  instrumentation  setup:  (a)  schematic  of  the  location  of  the  PWAS  sensors  on  the 

Aluminum  plate  as  well  as  the  equipment  necessary  to  generate  and  capture  the  signals;  (b) 
picture  of  the  actual  equipment  and  location  of  sensors 
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Table  3 


Location  of  PWAS  sensors  on  the  Aluminum  plate 


Sensor# 

X(mm) 

Y(mm) 

1 

100 

100 

2 

100 

175 

3 

100 

250 

4 

100 

325 

5 

100 

400 

6 

450 

100 

7 

450 

250 

8 

450 

400 

9 

800 

100 

10 

800 

250 

11 

800 

400 

The  excitation  signal  was  a  3-count  sine  burst  wave  with  Hanning  window.  The  Hanning  windowed 
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signal  is  defined  as:  =  0.5xdl  -  cos(<5;)] ,  where  co  =  —  ,  and  n  is  the  number  of  eounts  (Figure  38).  The 

n 

generated  signal  has  a  frequeney  of  300  kHz  (100  kHz  x  3  eounts)  for  the  So  mode  and  120  kHz  (40  kHz 
X  3  eounts)  for  the  Aq  mode. 


Figure  38  Excitation  signal 


5.3.2.  Results  and  Discussions 

The  results  are  presented  in  Figure  39  through  Figure  44,  whieh  show  the  transmitted  and  reeeived  signals 
using  pair  of  PWAS  sensors  in  a  reeiproeal  way  (#1  =  T,  #7  =  R;  #7  =  T,  #1  =  R);  (#8  =  T,#ll  =R;#ll  = 
T,  #8  =  R);  (#1  =  T,  #11  =  R,  #11  =  T,  #1  =  R).  First,  the  SO  mode  results  are  presented  (Figure  39 
through  Figure  41).  Then,  the  AO  mode  results  are  presented  (Figure  42  through  Figure  44).  Examining 
the  paired  signals  presented  in  Figure  39  through  Figure  44,  it  ean  be  eoneluded  that  reeiproeity  is 
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verified,  sinee  the  a)  and  b)  signals  in  every  figure  always  seem  identieal.  There  are  no  apparent 
differenees  in  the  time  of  arrival  or  in  the  amplitude  of  the  reeeived  signal  when  one  sensor  was  used  as 
transmitter  or  as  reeeiver.  This  proves  the  reeiproeity  relation  for  this  PWAS  installation  for  both  So  and 
Ao  modes  of  the  Lamb  wave. 


5.3.2. 1.  Results  for  the  So  mode 

PWAS  #1  PWAS  #7 


Figure  39  Transmitted  and  received  So  mode  Lamb  wave:  (a)  PWAS  #1  was  used  as  a  transmitter  and 
PWAS  #7  was  the  receiver;  (b)  PWAS  #7  was  used  as  a  transmitter  and  PWAS  #1  was  the 
receiver 


PWAS  #8  PWAS  #1 1 


Figure  40  Transmitted  and  received  So  mode  Lamb  wave:  (a)  PWAS  #8  was  used  as  a  transmitter  and 
PWAS  #1 1  was  the  receiver;  (b)  PWAS  #1 1  was  used  as  a  transmitter  and  PWAS  #8  was 
the  receiver 
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PWAS  #1  PWAS  #1 1 


Figure  41  Transmitted  and  received  So  mode  Lamb  wave:  (a)  PWAS  #1  was  used  as  a  transmitter  and 
PWAS  #1 1  was  the  receiver;  (b)  PWAS  #1 1  was  used  as  a  transmitter  and  PWAS  #1  was 
the  receiver 

5.3. 2.2.  Results  for  the  Aq  mode 


PWAS  #1  PWAS  #7 

Transmitter  Receiver 


Figure  42  Transmitted  and  received  Ao  mode  Lamb  wave:  (a)  PWAS  #1  was  used  as  a  transmitter  and 
PWAS  #7  was  the  receiver;  (b)  PWAS  #7  was  used  as  a  transmitter  and  PWAS  #1  was  the 
receiver 
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PWAS  #8  PWAS  #1 1 

Transmitter  Receiver 


Figure  43  Transmitted  and  received  Ao  mode  Lamb  wave:  (a)  PWAS  #8  was  used  as  a  transmitter  and 
PWAS  #1 1  was  the  receiver;  (b)  PWAS  #1 1  was  used  as  a  transmitter  and  PWAS  #8  was 
the  receiver 


Figure  44  Transmitted  and  received  Ao  mode  Lamb  wave:  (a)  PWAS  #1  was  used  as  a  transmitter  and 
PWAS  #1 1  was  the  receiver;  (b)  PWAS  #1 1  was  used  as  a  transmitter  and  PWAS  #1  was 
the  receiver 
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5.4.  Subtask  2.4  -  Modeling  and  Simulation  for  Detection  of  Corrosion  Using  Lamb 
Waves 

5.4.1.  Methods,  Assumptions,  and  Procedures 


Figure  45  Dispersion  characteristics  of  Aq  Lamb  wave  mode  in  1-mm  2024-T3  aluminum  plate:  (a)  Aq 
wave  speed  vs.  frequency;  (b)  Aq  wavelength  vs.  frequency 


Figure  46  Variation  of  tuned  excitation  frequency  with  PWAS  diameter  for  Aq  Lamb  wave  mode 

In  the  analysis,  the  exeitation  of  the  Aq  Lamb  wave  mode  whieh,  at  low  frequeneies,  resembles  flexural 
plate  waves.  The  highly  dispersive  eharaeteristie  of  Aq  flexural  waves  makes  them  very  sensitive  to  small 
ehanges  in  plate  thiekness,  as  expeeted  during  the  eorrosion  proeess.  Figure  45  shows  the  ealeulated 
dispersion  eharaeteristies  of  the  Aq  Lamb  mode  in  a  1-mm  2024-T3  aluminum  plate  in  the  1  to  100  kHz 
band.  It  is  notieed  that  the  wave  speed  and  wavelength  vary  rapidly  with  frequeney  in  this  low  frequeney 
band.  For  tuned  exeitation  with  a  10-mm  PWAS,  we  are  looking  for  wavelength  of  the  order  of  20  mm. 
Figure  45b  indieates  that  a  20-mm  wavelength  of  the  Aq  mode  would  be  eneountered  in  the  20-40  kHz 
range.  A  more  preeise  determination  of  the  PWAS  tuning  frequeney  is  given  in  Figure  46,  whieh  shows 
that  the  10-mm  diameter  PWAS  would  have  an  Aq  Lamb-mode  optimal  tuning  frequeney  of  around  30 
kHz.  Fine  tuning  experiments  deseribed  in  a  later  seetion  of  this  paper  have  found  the  exaet  tuning 
frequeney  of  31.5  KHz.  The  differenee  between  the  theoretieal  predietion  and  experimental  results  is 
attributable  to  the  2-D  effeets,  whieh  were  not  eaptured  in  this  1-D  analysis. 


38 


The  dispersive  propagation  of  the  Aq  Lamb- wave  mode  was  analyzed  using  the  Fourier  transform 
method.  The  exeitation  eonsisted  of  a  3.5-eount  Hanning  windowed  smoothed  tone  burst  of  frequeney  fo  = 
31.5  kHz  starting  at  time  to  and  lasting  for  the  duration  T: 


p(t)= 


,  if  tf)  ^  t  tf)  T 


0  otherwise 


(32) 


where  cOq  The  tone  burst  duration  was  ealeulated  in  terms  of  the  number  of  eounts,  Ncounts  =  3.5 

using  the  expression: 


T  = 


fo 


(33) 


The  exeitation  was  transformed  to  the  frequeney  domain  using  the  Fourier  transform: 


p{f)=rpit)e-^^^dt 

J— 00 


(34) 


In  the  frequeney  domain,  the  exeitation  was  multiplied  with  the  wave  transfer  funetion  (Doyle,  1997): 


(35) 


to  generate  the  frequeney  domain  response: 

u{co,x)  =  p{co)  •  G{co,x) 


(36) 


The  frequeney-dependent  wave  number  k{co)  of  Equation  (35)  was  ealeulated  using  the  low- frequeney 
approximation  of  the  Aq  wave  speed,  i.e.. 


k{(D)  = 


3p[l-v^ 

d^E 


CD 


1/2 


(37) 


Applying  the  inverse  Fourier  transform  to  Equation  (36),  we  obtained  the  Aq  mode  Lamb  wave  response 
as  a  funetion  of  time,  t,  and  loeation  in  the  plate,  x,  i.e.. 


I  P  +00 

(t,x)  =  —  u{o),x)e  da> 


(38) 


5.4.2.  Results  and  Discussions 

Figure  47  presents  the  simulation  of  the  Aq  mode  Lamb  wave  response  at  the  PWAS  souree,  x  =  0,  and  at 
distanees  x  =  250V2  and  x  =  500V2  mm  from  the  PWAS  souree.  It  is  apparent  that,  as  the  distanee  from 
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the  source  increases,  the  dispersion  of  the  Aq  mode  Lamb  wave  packet  increases.  Figure  48  presents  the 
effect  of  plate  thickness  variation  on  the  wave  packet  dispersion.  Under  consideration  is  a  wave  packet 
starting  at  tq  =  0  and  arriving  at  location  ri  =  250V2  (Figure  48a)  while  traveling  through  pristine  plate,  or 
at  ri.cor  =  250V2  (Figure  48b)  while  traveling  through  corroded  plate  that  has  experienced,  on  average,  a 
thickness  reduction  of  25%.  Comparison  of  Figure  48a  with  Figure  48b  reveals  that  the  effect  of  corrosion 
slows  the  wave  packet  down  and  to  increase  its  dispersion.  We  notice  that  the  packet  traveling  through  the 
pristine  region  seems  to  reach  its  peak  amplitude  at  around  280  micro-sec,  while  the  packet  traveling 
through  the  region  with  25%  average  corrosion  seems  to  reach  its  peak  amplitude  at  around  320  micro¬ 
sec.  This  40  micro-sec  delay  represents  an  average  increase  of  13%  in  travel  time.  Results  similar  to  these 
predictions  are  to  be  found  in  the  experimental  section  of  this  paper. 


Figure  47  Simulation  of  Aq  mode  Lamb  wave  response  at  x  =  250V2  and  500V2  mm  from  the  PWAS 
source  placed  at  x  =  0  mm  (3.5-count  smoothed  tone  burst  of  31 .5  kHz,  plate  thickness  h  =  1 
mm) 


Figure  48  Simulation  of  Aq  Lamb  wave  mode  response  at  250V2  mm  from  the  PWAS  source  for  two 
plate  thickness  values:  (a)  h  =  1  mm;  (b)  =  0.75  mm  (3.5-count  smoothed  tone  burst  of 

31.5  kHz) 

Figure  49  presents  signal  cross  plots  between  the  signal  simulated  for  a  pristine  plate  (“pristine  signal”) 
and  the  signal  simulated  for  a  corroded  plate  (“corroded  signal”).  These  plots  are  very  sensitive  to  small 
shifts  in  signal  phase.  Two  charts  are  presented,  one  for  0.5%  total  material  loss  (Figure  49a)  and  one  for 
1.0%  total  material  loss  (Figure  49b).  In  each  chart,  two  plots  are  made:  first,  a  plot  of  the  pristine  signal 
on  itself,  which  results  in  a  straight  45®  line;  second,  a  plot  of  the  corroded  signal  against  the  pristine 
signal.  This  second  plot  results  in  an  ellipse.  By  comparing  Figure  49a  with  Figure  49b  we  observe  that 
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the  width  of  the  ellipse  eorrelates  with  the  simulated  eorrosion  intensity.  For  mild  eorrosion  (0.5%  total 
material  loss),  the  ellipse  is  relatively  slim  (Figure  49a).  For  more  intense  eorrosion,  (1.0%  total  material 
loss)  the  ellipse  is  mueh  wider.  These  results  indieate  that  the  signal  eross  plots  ean  be  effeetively  used  to 
identify  the  signal  phase  shifts  and  dispersions  due  to  the  presenee  of  eorrosion. 


0.5%  total  material  loss 


Pristine  signal,  au 


1.0%  total  material  loss 


Figure  49  Cross  plot  between  the  pristine  and  corroded  signals:  (a)  0.5%  total  material  loss;  (b)  1% 
total  material  loss 


5.5.  Task!  -  Conclusions 

The  purpose  of  Task  2  was  the  find  the  optimum  design  for  the  PWAS.  PWAS  ean  be  used  in  many  ways 
for  struetural  health  monitoring  and  non-destruetive  evaluation,  but  if  the  design  of  the  PWAS  is  not 
optimized,  the  effeetiveness  and  performanee  of  the  PWAS  ean  be  eompromised. 
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In  this  task  we  started  with  the  study  of  the  Lamb  waveforms  and  surfaee  displaeements  for  various 
material  properties  and  geometries  (thiekness,  size,  ete).  Theoretieal  predietions  of  the  PWAS’  interaetion 
with  Lamb  waves  were  used  as  a  starting  base  for  optimizing  the  design  and  performanee  of  future  health 
monitoring  systems.  We  identified  the  effeet  of  frequeney  and  sensor  size/shape  on  the  sensor 
effeetiveness.  These  are  important  parameters  that  eontribute  to  the  design  of  a  good  health  monitoring 
system. 

Mode  tuning  software  simulation  was  developed  to  be  able  to  seleetively  tune  into  speeifie  Lamb  modes. 
It  has  been  observed  in  separate  studies  that  the  deteetion  of  through-the-thiekness  eraeks  with  the  pulse- 
eeho  method  is  mueh  better  with  the  So  mode  than  with  the  Aq  mode.  On  the  other  hand,  the  Aq  mode 
seems  to  better  suited  for  the  deteetion  of  delaminations,  disbonds,  and  eorrosion  with  the  piteh-eateh 
aeousto-ultrasonies  teehniques.  The  mode-tuning  software  developed  in  this  task  faeilitates  the  optimum 
design  of  PWAS  transdueers  to  aehieve  the  desired  damage  deteetion  results. 

Modeling  using  MATLAB  and  MathCAD  was  used  to  simulate  the  Lamb  waves  and  prediet  their 
behavior.  This  modeling  was  used  to  find  the  optimum  PWAS  design  for  various  applieations.  The  sensor 
design  eriteria  were  tested  experimentally.  In  our  experiments  we  verified  the  reeiproeity  of  PWAS  Lamb 
wave  transmission  and  reeeption  and  determined  that  optimum  tuning  into  Aq  and  So  Lamb  wave  modes 
ean  be  aehieved.  The  signals  predieted  by  the  theoretieal  modeling  behaved  similar  to  the  signals 
measured  in  the  aetual  test.  This  ean  be  used  to  optimize  the  sensor  design  by  understanding  how  the 
waves  behave  during  their  travel  in  the  plate  material.  Comparative  presentation  of  predieted  and 
measured  signals  indieated  that  a  reliable  predietion  method  has  been  obtained  that  ean  be  used  for 
optimum  sensor  design. 
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6.  TASK  3:  STEERING-BEAM  LARGE  AREA  COVERAGE 

This  task  will  extend  the  preliminary  work  done  on  embedded  ultrasonie  struetural  radar  (EUSR).  To 
date,  the  EUSR  eoneept  has  been  used  to  demonstrate  the  use  of  the  seanning  beam  to  deteet: 

■  20  mm  horizontal  eraek  at  90°  (broadside  loeation,  i.e.  vertieal  above  the  array) 

■  20  mm  horizontal  eraek  at  135°  (offside  loeation,  i.e.  to  the  left  of  the  array) 

Both  eraeks  were  positioned  at  a  radial  distanee  of  300  mm  from  the  array.  Thin-plate  speeimens  1220- 
mm  sq.,  1-mm  thiek  2024  T3  were  used.  In  the  new  tests,  we  plan  to  aehieve  the  following:  examine 
multiple  eraeks;  the  effeet  of  eraek  orientation;  the  effeets  of  eraek  size  and  imaging  large  eraeks;  and  the 
effeets  of  eurvature  -  theoretieally,  guided  waves  are  able  to  propagate  in  a  plate  with  small  eurvature  and 
verify  this  for  various  eurvatures  and  identify  the  severe  eases. 


Task  Overview 

■  Subtask  3.1-  EUSR  With  Curvature 

■  Subtask  3.2  -  The  Effeets  of  Craek  Orientation  and  Multiple  Targets  on  EUSR 

■  Subtask  3.3  -  Automatie  Signal  Aequisition 


6.1.  Subtask 3.1  -EUSR  With  Curvature 

6.1.1.  Methods,  Assumptions,  and  Procedures 

Experiments  were  done  on  speeimens  with  three  different  extents  eurvature.  To  faeilitate  the  experiments, 
the  extent  of  eurvature  was  defined  by  L,  the  differenee  between  the  original  length  L  and  the  ehord  C, 
See  Figure  50  below.  We  went  through  5mm,  10mm,  15mm  and  20mm  to  eover  the  three  extents. 
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Figure  50  -  Curvature  of  the  speeimen 


During  the  proeess  to  set  up  the  eurvature,  it  was  diseovered  that  there  eould  be  two  types  of  eurvature  on 
the  square  plate.  One  is  parallel  to  the  PWAS  array;  the  other  is  perpendieular  to  the  PWAS  array.  The 
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parallel  PWAS  array  will  be  ealled  ‘direetion  1’  and  the  perpendieular  PWAS  array  will  be  ealled 
‘direetion  2’  (Figure  51). 


6.1.2.  Results  and  Discussions 

In  the  experiments,  both  direetions  were  tested  to  verify  that  guided  waves  are  able  to  be  used  as  the 
propagating  wave  sinee  its  propagation  is  eonfined  within  the  strueture  or  material.  Therefore  the 
eurvature  should  not  affeet  the  deteetability  of  embedded  ultrasonie  struetural  radar  (EUSR). 

Aeeording  to  the  experimental  results,  it  was  eoneluded  that  generally  the  eurvature  will  not  affeet  the 
effieieney  of  EUSR  deteetability. 


Figure  51  -  Two  types  of  curvature:  a)  direction  1,  parallel  to  the  PWAS  array;  b)  direction  2, 
perpendicular  to  the  PWAS  array. 
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6.1. 2.1.  Experiment  with  0  curvature 


Figure  52  -  Experiment  with  0  curvature 


Figure  53  -  EUSR  analysis  for  0  curvature 
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6.1. 2.2.  Experiment  with  5mm  shortening  (Curvature  radius  R  =  3900  mm)  in  direction  1 


Figure  54  -  Experiment  with  5mm  curvature  at  direction  1 


EUSR_V3.Vi 


Figure  55  -  EUSR  analysis  for  5mm  curvature  at  direction  1 
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6.1. 2.3.  Experiment  with  10mm  shortening  (curvature  radius  R  =2750  mm)  in  direction  1 


Figure  56  -  Experiment  with  10mm  curvature  at  direction  1 


Figure  57  -  EUSR  analysis  for  10mm  curvature  at  direction  1 
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6.1. 2.4.  Experiment  with  15mm  shortening  (curvature  radius  R  =  2240  mm)  in  direction  1 


Figure  58  -  Experiment  with  15mm  curvature  at  direction  1 


^  EU5R_V3.Vl 


Figure  59  -  EUSR  analysis  for  15mm  curvature  at  direction  1 
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6.1. 2.5.  Experiment  with  20mm  shortening  (curvature  radius  R  =  1950  mm)  in  direction  1 


Figure  60  -  Experiment  with  20mm  curvature  at  direction  1 


Figure  61  -  EUSR  analysis  for  20mm  curvature  at  direction  1 
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6.1. 2.6.  Experiment  with  5mm  shortening  (curvature  radius  R  =  3900  mm)  in  direction  2 


Figure  62  -  Experiment  with  5mm  curvature  at  direction  2 


^  EUSR_v3,\Ji 


Figure  63  -  EUSR  analysis  for  5mm  curvature  at  direction  2 
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6.1.2. 7.  Experiment  with  10mm  shortening  (curvature  radius  R  =2750  mm)  in  direction  2 


Figure  64  -  Experiment  with  10mm  curvature  at  direction  2 


W,  EU5R_v3,vi 


Figure  65  -  EUSR  anaiysis  for  10mm  curvature  at  direction  2 
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6.1. 2.8.  Experiment  with  15mm  shortening  (curvature  radius  R  =  2240  mm)  in  direction  2 


Figure  66  -  Experiment  with  15mm  curvature  at  direction  2 


EU5R_v3,vi 


Figure  67  -  EUSR  analysis  for  15mm  curvature  at  direction  2 
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6.1. 2.9.  Experiment  with  20mm  shortening  (curvature  radius  R  =  1950mm)  in  direction  2 


Figure  68  -  Experiment  with  20mm  curvature  at  direction  2 


EUSR_u3.vi 


Figure  69  -  EUSR  analysis  for  20mm  curvature  at  direction  2 
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6.2. 


Subtask  3.2  -  The  Effects  of  Crack  Orientation  And  Multiple  Targets  onEUSR 


6.2.1.  Methods,  Assumptions,  and  Procedures 

6.2.1. 1.  The  Embedded  Ultrasonic  Structural  Radar  (EUSR)  Algorithm  and  Previous  Work 

The  principle  of  operation  of  the  embedded  ultrasonic  structural  radar  (EUSR)  is  derived  from  two 
general  principles:  (1)  The  principle  of  guided  Lamb  wave  generation  with  piezoelectric  wafer  active 
sensors  (PWAS);  (2)  The  principles  of  conventional  phased-array  radar. 

The  guided  Lamb  waves  generated  by  PWAS  have  the  important  property  that  they  stay  confined  inside 
the  walls  of  a  thin-wall  structure,  and  hence  can  travel  over  large  distances.  In  addition,  the  guided  wave 
can  also  travel  inside  curved  walls,  which  makes  them  ideal  for  applications  in  the  ultrasonic  inspection 
of  aircraft,  missiles,  pressure  vessel,  oil  tanks,  pipelines,  etc.  This  was  proved  by  curvature  experiment. 
Lamb  waves  can  exist  in  a  number  of  dispersive  modes.  However,  through  smoothed  tone -burst  excitation 
and  frequency  tuning,  it  is  possible  to  confine  the  excitation  to  a  particular  Lamb  wave  mode,  of  carrier 
frequency  fc,  wave  speed  c,  and  wave  length  X  =  cl  .  Hence,  the  smoothed  tone-burst  signal  generated 
by  one  PWAS  is  of  the  form: 

Sj,  {t)  =  {t)  cos  It:  f^t  (39) 

where  ^o(0  is  a  short-duration  smoothing  window  that  is  applied  to  the  carrier  signal  of  frequency  fc 
between  0  and  tp. 

The  principles  of  conventional  phased-array  radar  are  applied  to  the  PWAS-generated  guided  waves, 
assuming  a  uniform  linear  array  of  M  sensors  (PWAS),  with  each  PWAS  acting  as  a  point  wise  omni¬ 
directional  transmitter  and  receiver.  The  PWAS  in  the  array  are  spaced  at  the  distance  d,  which  is 
assumed  much  smaller  than  the  distance  r  to  a  generic,  far-distance  point,  P.  Since  d«r,  the  rays  joining 
the  sensors  with  the  point  P  can  be  assimilated  with  a  parallel  fascicle,  of  angle  (j)  .  Therefore,  for  the 

PWAS,  the  distance  will  be  shorted  by  m{d  cos^?^)  .  If  all  the  PWAS  are  fired  simultaneously,  the  signal 
from  the  mth  PWAS  will  arrive  at  P  quicker  by  A^(^)  =  m(d  cos^?))/  c  .  Yet,  if  the  PWAS  are  not  fired 
simultaneously,  but  with  some  individual  delays,  m  =  0,1,..., M-1,  then  the  total  signal  received  at 
point  P  will  be: 

I  ^-1  (  r  \ 

=  -p  Z  t - +  (40) 

Yr  ^=0  \  c  J 

Xjr  represents  the  decrease  in  the  wave  amplitude  due  to  the  omni-directional  2-D  radiation,  and  rjc  is 
the  delay  due  to  the  travel  distance  between  the  reference  PWAS  {m  =  0)  and  the  point  P.  (Here  wave- 
energy  conservation,  i.e.,  no  dissipation,  is  assumed.) 

Transmitter  beam  forming:  if  we  have  5^  =  mA(^)  ,  then  Equation  (40)  becomes: 


s^{t)=  M  -  t  — 


That’s  to  say,  there  is  an  M  times  increase  in  the  signal  strength  with  respect  to  a  simple  sensor.  This 
leads  directly  to  the  beam  forming  principle  that  if  5^  =  md  cos((^q)/  c  ,  and  since  A^  =  md  cos((^)/ c  , 

then  constructive  interference  (beam  forming)  takes  place  when  cos(^2J)  =  cos(^2Jq)  ,  i.e.  at  angles  ^  = 
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and  (j)  =  .  Thus,  the  forming  of  a  beam  at  angles  and  is  aehieved  through  delays  in  the  firing 

of  the  sensors  in  the  array. 


Reeeiver  beam  forming:  if  the  point  P  is  an  omni-direetional  souree  at  azimuth  then  the  signals 
reeeived  at  the  m^^  sensor  will  arrive  quieker  by  m^cf)  =  md  cos(^q)/  c  .  Henee,  we  ean  synehronize  the 
signals  reeeived  at  all  the  sensors  by  delaying  them  by: 


5 

m 


(<^o)  =  m-cos(<zio) 


(42) 


Pulse-eeho  method:  Assume  that  a  target  exists  at  azimuth  and  distanee  R.  The  transmitter  beam 
former  is  sweeping  the  range  in  inereasing  angles  (f)  and  reeeives  an  eeho  when  (/)  =  (j)^.  The  eeho  will  be 

reeeived  on  all  sensors,  but  the  signals  will  not  be  synehronized.  To  synehronize  the  sensors  signals,  the 
delays  5^  =  md  cos((2Jq)  /  c  need  to  be  applied. 


The  original  signal  is: 


5p(0  = 


M 

7r‘ 


2R^ 

t - 

c  j 


(43) 


At  the  target,  the  signal  is  backscattered  with  a  backscatter  coefficient,  A.  Hence,  the  signal  received  at 
each  sensor  will  be: 


A-M 

R 


f 


2R 


\ 


^ - 1"  ((^) 

V  c 


(44) 


The  receiver  beam  former  assembles  the  signals  from  all  the  sensors  with  the  appropriate  time  delays: 


A-M 

^r(0  =  — 

K  m=0 


2  R 

t - + 


(45) 


Constructive  interference  between  the  received  signals  is  aehieved  when^J^  =  md  cos(^q)/  c  .  Thus,  the 
assembled  reeeive  signal  will  be  again  boosted  M  times,  with  respeet  to  the  individual  sensors: 


A-M 

^r(0  = - 

R  “=« 


2R'] 

t - 

I  c  J 


(46) 
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Figure  70  the  transmitted  signal  and  the  received  signal 
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In  general,  a  target  eraek  is  unknown,  i.e.  the  target  loeation  is  unknown.  Sinee  we  use  the  polar 
eoordinates  in  the  radar  system,  the  loeation  of  an  unknown  target  is  defined  by  the  angle  (j)^ .  The  eoarse 

estimation  of  (j)^  is  implemented  by  using  the  (j)^  sweeping  method.  EUSR  will  sean  through  0°  to  180° 
by  inerementing  (j)^  by  1°  eaeh  time,  until  the  maximum  reeeived  energy  is  obtained.  max£’j^(^Q)  is  the 
maximum  reeeived  energy  by  the  definition: 

=  1^^'“  I  P  dt  (47) 

P 

Figure  70  is  an  example  of  the  transmitted  signal  and  reeeived  signal  pair  eaptured  in  our  experiment. 

The  speeimens  we  used  in  the  lab  are  4-feet  long,  1-mm  thiek  square  aluminum  plates  with  an  8-element 
PWAS  array  at  the  eenter.  The  layout  is  shown  in  Figure  71  (a).  Figure  71  (b)  is  the  experiment  setup. 
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(24.00  in) 
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d  =  1  mm  8  round  PWAS  array, 
1mm  apart  , 


(a) 


1220-mm  sq.,  1-mm  thick  2024  T3 
(22.45-in  sq.,  0.040-in  thick) 


Figure  71  EUSR  experiment  with  the  8  element  PWAS  array;  (a)  specimen  layout;  (b)  experiment  setup 


In  the  experiment,  broadside  position  is  defined  to  be  perpendieular  to  the  PWAS  array  and  offside 
position  is  anywhere  between  the  0°  and  180°  but  the  broadside  position.  The  eraeks  are  simulated  by 
through-plate  eraeks  14mm  long  and  0.5mm  wide. 

Previous  proof-of-eoneept  work  has  been  done  by  using  FUSR  to  verify  its  ability  for  deteeting  a  single 
eraek  at  either  broadside  or  offside  position  on  a  thin-plate  speeimen.  Figure  72  (b)  and  Figure  73  (b)  are 
the  2D  display  of  single  eraek  inspeetion  by  using  our  previous  FUSR.  The  offside  position  is  angle  137°. 


Figure  72  single  broadside  crack  detection  with  previous  EUSR:  (a)  schematic;  (b)  EUSR  GUI  mapped 
image 
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o.e 


(a) 

Figure  73  single  offside  crack  detection  with  previous  EUSR:  (a)  schematic;  (b)  EUSR  GUI  mapped 
image 

6,2, 1,2,  An  Improved  EUSR  by  using  the  Hilbert  Transform 

From  the  specimen  inspection  plots,  it  was  determined  that  they  can  generally  locate  the  damage;  they  are 
far  than  close  to  the  real  cracks.  A  solution  to  improve  the  image  quality  and  precision  was  needed  to  be 
developed.  It  was  noticed  that  echoes  from  the  damage  do  not  have  single  peaks  due  to  the  dispersion 
property  or  some  reasons.  Something  had  to  be  done  to  get  the  echoes  close  to  the  theoretical  ones,  having 
single  peaks.  The  method  that  was  explored  was  the  Hilbert  transform,  which  is  able  to  extract  the 
envelope  of  a  curve. 

The  envelope  is  a  curve  or  surface  that  is  tangent  to  every  one  of  a  family  of  curves  or  surfaces.  Envelope 
extracts  the  amplitude  of  a  periodic  signal.  It  can  be  used  to  simplify  the  process  of  detecting  the  time  of 
arrival  for  the  wave  packets  in  our  EUSR  system.  In  the  EUSR,  the  envelope  of  the  signal  is  extracted  by 
applying  Hilbert  transform  to  the  cross  correlation  signal. 

The  Hilbert  transform  is  defined  as: 


0.f5 

0.3 

0.2 

0.1 


(b) 


...  , , ,  1  x(t)  , 

= - }  -^dr  (48) 

7T  — 00  I  ^ 

Hilbert  transform  can  also  be  used  to  construct  a  complex  signal: 


V0  =  %e(0  +  ^-^Im(0 


Where 


%e(0  =  ^(0 

^Im(O  =  ^(^(0) 


(49) 


(50) 
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Figure  74  a  tone-burst  signal  and  its  envelope  extracted  by  using  Hilbert  transform 


(a)  me  ter  (b)  meter 

Figure  75  -  single  crack  detection  with  the  improved  EUSR:  (a)  broadside  crack  detection;  (b)  offside 
crack  detection 

The  real  part  of  the  eonstrueted  signal,  ?  is  the  original  data  x{t) ,  while  the  imaginary  part 

is  the  Hilbert  transform  of  x{t) .  Aetually,  the  imaginary  part  is  a  version  of  the  original  real  sequenee 
after  a  90°  phase  shift.  Thus,  the  Hilbert  transformed  signal  has  the  same  amplitude  and  frequeney  eontent 
as  the  original  real  signal  and  ineludes  phase  information  that  depends  on  the  phase  of  the  original  signal. 
The  magnitude  of  eaeh  eomplex  value  has  same  the  amplitude  as  the  origin  signal.  Therefore,  the 
magnitude  of  the  analytieal  signal  is  the  envelope  of  the  original  signal.  Just  by  observing  the  envelope 
signal,  the  wave  paekages  ean  be  easily  reeognized.  Figure  74  is  an  example  showing  the  envelope 
extraeted  after  applying  the  Hilbert  transform. 

The  proeessing  results  show  that  this  proeess  ean  result  in  a  more  distinet  representation  of  eehoes  from 
the  eraek  and  therefore  a  more  elear  indieation  of  the  eraek.  Compared  with  Figure  73,  Figure  75  shows 
results  of  the  inspeetion  of  the  eraek  deteetion  with  better  resolution  and  more  elose  to  the  eraek  real  size. 
Meanwhile,  the  new  eraeks  inspeetion  has  fewer  disturbanees  from  the  side-lobe  effeets. 

Therefore,  the  Hilbert  transform  will  be  apply  after  the  EUSR  algorithm  in  order  to  extraet  the  envelopes 
of  the  reeonstrueted  EUSR  signals  sueh  that  an  eeho  ean  be  ereated  with  a  single  peak.  Figure  76  shows 
the  data  flow  of  the  EUSR. 
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Figure  76  -  data  flow  of  the  damage  detection  by  using  the  improved  EUSR  algorithm 
6.2.1. 3.  EUSR  multi-crack  detection  experiments 

The  previous  work  has  proved  EUSR’s  ability  to  deteet  a  single  through  plate  eraek  in  a  thin  plate.  The 
new  experiments  are  designed  to  verify  EUSR  is  also  able  to  deteet  multiple  eraeks.  The  following 
experiments  were  preformed: 

•  Broadside  crack  with  certain  slope 

•  Two  aligned  offside  cracks 

•  Three  aligned  cracks,  two  offside  cracks  and  one  horizontal  broadside  crack 

•  Three  aligned  cracks,  two  offside  cracks  and  one  perpendicular  broadside  crack 

6.2.2.  Results  and  Discussions 


6. 2. 2.1.  Broadside  Crack  with  Certain  Slope 


The  speeimen  #7  has  a  broadside  eraek  with  the  20-mm  eraek  dimension,  but  30°  to  the  horizontal.  The 
experiment  result  in  Figure  77  shows  that  EUSR  is  able  to  indieate  the  slope.  The  plot  is  different  from 
the  horizontal  eraek  showing  in  Figure  75  a).  In  the  future,  further  experimentation  will  be  done  to  find 
out  the  minimum  slope  that  EUSR  is  sensitive  and  able  to  indieate,  and  eompare  the  results  of  different 
slopes. 
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6,2,2,2,  Two  aligned  offside  cracks 

The  specimen  #6  has  two  symmetrical  offside  cracks  at  63®  and  117®  positions,  respectively.  Figure  78 
shows  the  EUSR  inspection  result. 

One  thing  has  to  be  mentioned  is,  for  some  reason,  the  two  cracks  are  actually  not  exactly  symmetrical  to 
each  other  as  the  original  design.  They  are  3  mm  from  each  vertically.  However,  the  2D  inspection  plot 
picked  out  such  a  difference,  and  it  is  shown  that  the  two  cracks  are  in  different  circles  (having  different 
vertical  distance  to  the  PWAS  array).  It  can  be  concluded  that  the  EUSR  is  able  to  even  pick  up  that  small 
difference  (Figure  78). 


0.6 

0.6 

0.3 

0.2 

0.1 


V 

> 

(b) 


meter 


Figure  78  -  the  two  symmetrical  offside  cracks  at  63°  and  137°,  respectively:  (a)  schematic;  (b)  EUSR 
GUI  mapped  image 


6,2,23.  Next  Step:  Three  Aligned  Cracks 

In  the  three  aligned  cracks  experiments,  two  designs  were  created.  Both  designs  have  two  offside 
horizontal  cracks  at  63®  and  137®  positions.  Design  1  has  a  broadside  crack  parallel  to  the  PWAS  array 
while  design  2  has  a  broadside  crack  perpendicular  to  the  PWAS  array  (Design  1:  the  horizontal 
broadside  crack  and  Design  2:  the  vertical  broadside  crack). 


Figure  79  -  three  aligned  cracks:  (a)  schematic  of  the  horizontal  broadside  crack;  (b)  schematic  of  the 
vertical  broadside  crack 

To  this  moment,  EUSR  algorithm  is  able  to  detect  the  presence  of  the  three  cracks  in  the  A-scan  signal, 
but  not  in  the  mapped  EUSR  image  for  the  three  cracks  experiments.  The  reason  for  this  may  be  that  the 
broadside  crack  reflects  more  energy  than  the  offside  cracks  and  therefore  the  energy  reflected  from  the 
perpendicular  crack  dominates  the  EUSR  image.  For  the  design  2,  recalling  the  result  attained  from  the 
pin-hole  experiments  (Task  4),  the  minimum  hole  size  that  can  be  detected  by  EUSR  is  1.57  mm. 
Therefore  the  crack  in  the  design  2  may  be  too  tiny  to  be  detected. 
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In  the  next  stage,  the  signal  proeessing  teehniques  will  be  applied  and  signal  analysis  methods  will  be 
used  to  improve  the  results.  The  vertieal  eraek  will  gradually  be  enlarged  to  find  out  the  smallest 
dimension  deteeted  and  further  verify  EUSR. 


6.3.  Subtask  3.3 -Automatic  Signal  Acquisition 

6.3.1.  Methods,  Assumptions,  and  Procedures 

An  automatie  signal  switch  is  necessary  to  send  the  excitation  signals  to  PWAS  and  acquire  the  response 
signal  from  another  PWAS.  Such  a  program-controlled  switch  can  quickly  and  precisely  execute  the  data 
collection  in  a  way  which  is  more  efficient  and  reliable  than  the  manual  switching  operations.  The 
automation  of  data  collection  consists  of  two  parts:  the  hardware  part  of  auto  signal  switch  box  and  the 
corresponding  software  part  of  PC  control  program.  In  this  method,  digit  control  signals  are  generated  by 
the  PC  software  and  sent  to  the  switch  box  through  the  parallel  port.  According  to  the  control  signals  form 
the  parallel  port,  the  switch  box  will  connect  the  function  generator  and  oscilloscope  each  to  one  sensor 
(these  two  sensor  can  be  the  same)  of  the  PWAS  array  respectively.  Thus,  one  signal  measurement  route 
is  constructed,  the  excitation  signal  is  transmitted  to  the  PWAS  array  and  echo  signals  are  received  by  the 
oscilloscope.  With  this  method,  the  measurement  loops  are  performed  automatically  under  the  control  of 
the  PC  software. 


t' 


Auto  signal  switch  unit 

Figure  80  Experimental  setup 

Figure  80  shows  the  whole  measurement  equipments  setup.  The  PWAS  array  is  connected  to  the  switch 
box  with  an  8-pin  ribbon  bus;  the  function  generator  and  oscilloscope  are  connected  to  the  switch  box 
with  coaxial  cables.  The  switch  box  is  connected  to  the  parallel  port  of  the  control  PC  to  receive  digit 
control  signals. 
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The  software  part,  developed  in  Lab  VIEW,  eontrols  the  hardware  part.  The  “out  port”  fiinetion  in 
Lab  VIEW  is  used  to  send  digit  signals  through  PC  parallel  port.  This  software  provides  a  graphie  user 
interfaee  (GUI)  to  faeilitate  the  data  eolleetion  of  PWAS  array  (Figure  81).  With  the  GUI,  user  ean 
eonfigure  the  switeh  unit  work  in  an  auto  signal  aequiring  mode  in  whieh  signal  transmitted  to  and 
reeeived  from  assigned  sensors  ean  be  eompleted  automatieally  without  ehanging  the  hardware 
eonneetion  by  hand.  When  the  auto  switeh  is  in  the  automatie  mode,  user  will  input  2  numbers  and  the 
path  name,  the  auto  switeh  will  perform  the  measurement  loops  that  start  from  the  first  number  until  the 
seeond  number  and  the  data  from  these  measurement  loops  will  be  saved  in  that  path;  when  in  manual 
mode,  the  auto  switeh  will  allow  the  user  to  eolleet  data  with  the  transmitting  and  reeeiving  sensors 
speeified  by  the  two  number  inputs.  After  these  parameters  defined,  the  eontrol  software  will  send  out  8- 
bit  digit  signals  through  parallel  port  and  these  signals  will  then  be  deeoded  to  eontrol  the  read-relays. 

Two  rows  of  indieating  LEDs  will  be  lit  in  green  eolors  to  show  whieh  sensor  is  transmitting  exeitation 
signals  and  whieh  one  is  used  to  reeeive  eeho  signals.  During  the  data  eolleetion  proeess,  the  wave  form 
will  also  be  displayed  on  this  GUI. 


The  eontrol  program  is  easy  to  implement  and  ean  be  integrated  into  an  upper  level  program  that  exeeutes 
the  whole  task  of  signal  aequisition  and  analysis.  Beeause  of  the  eoneise  design  of  the  hardware,  the 
eoneept  of  the  auto  signal  switeh  ean  be  extended  to  other  applieation  eases  sueh  as  the  eleetromeehanieal 
(E/M)  impedanee  measurement  for  SHM. 


Figure  81  Autoswitch  device  for  automated  signal  acquisition  in  the  steering  beam  experiments:  (a) 
experimental  setup;  (b)  graphical  user  interface  of  the  controlling  LabVIEW  program 


6.4.  Task  3-  Conclusions 

In  Subtask  3.1,  according  to  the  experimental  results,  it  was  concluded  that  generally  the  curvature  will 
not  affect  the  efficiency  of  EUSR  detectability.  Subtask  3.2  has  presented  the  applications  of  signal 
processing  techniques  to  the  existing  embedded  ultrasonic  structural  radar  (EUSR)  algorithm  using 
piezoelectric  wafer  active  sensors.  After  presenting  the  general  principles  of  the  EUSR  algorithm  and  the 
previous  work,  a  better  method  was  presented  (the  Hilbert  transform)  to  extract  the  signal  envelopes  to 
improve  the  efficiency  and  reliability  of  existing  EUSR.  A  series  of  experiments  are  conducted  to  verify 
the  EUSR’s  ability  to  detect  multi-cracks,  including:  (a)  broadside  crack  with  certain  slope;  (b)  two 
horizontal  aligned  cracks.  The  inspection  plots  of  the  cracks  from  the  specimens  were  achieved  under 
monitoring.  For  the  three  horizontal  aligned  cracks  with  two  types  of  middle  crack  cases,  it  is  envisioned 
to  apply  signal  processing  techniques  and  enlarge  the  vertical  crack  dimension  in  our  next  step  coming 
soon.  In  Subtask  3.3,  the  automatic  signal  switch  provided  fast  data  acquisition  from  arrays  and  a 
simple  solution  for  sensor  array  connection  and  data  collection. 
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7.  TASK  4:  PROBABILITY  OF  DETECTION  (POD)  WITH  PWAS  PROBES  (AN  INITIAL 
ANALYSIS) 

The  aging  process  of  commercial  and  military  fleet  has  generated  an  increasing  attention  to 
nondestructive  testing  (NDT)  methods.  Because  problems  in  aging  aircrafts  are  growing  faster  than  the 
rate  at  which  aircraft  age,  it  is  important  to  transfer  these  NDT  techniques  from  the  laboratory  to  the 
aircraft  as  soon  as  possible.  Part  of  this  transfer  is  the  probability  of  detection  (POD)  of  anomalies.  Before 
implementing  NDT  methods  on  the  aircraft  it  is  important  to  know  if  the  methods  can  detect  anomalies 
and  what  would  be  the  efficiency  of  detection. 

One  definition  of  the  POD  as  a  function  of  the  damage  size  is  the  fraction  of  damages  of  a  nominal  size 
that  are  expected  to  be  detected.  This  can  be  expressed  as  a  ratio  of  the  probability  of  detecting  damage 
with  a  certain  confidence  interval  (usually  90%,  95%).  The  information  obtained  in  this  way  is  usually 
presented  in  a  graphical  form,  where  the  probability  of  detection  is  plotted  as  a  function  of  the  damage  for 
a  fixed  confidence  interval  as  presented  in  Figure  82  (Grills,  2001). 


Figure  82  Typical  POD  Curves  (Grills,  2001) 

It  is  also  important  to  keep  the  probability  of  false  alarm  (POFM)  below  a  certain  level.  A  false  alarm 
occurs  when  an  anomaly  is  detected  when  in  fact  there  is  no  anomaly.  In  real  applications  it  is  desired  not 
to  miss  any  damages  so,  most  detection  algorithms  are  set  to  insure  that  damage  is  not  missed  while 
keeping  the  false  rate  below  1%. 

7.7.  Methods,  Assumptions,  and  Procedures 

7.1.1.  Planning  the  POD  Experiment 

The  purpose  of  this  work  is  to  identify  damage  in  simulated  aircraft  panels  using  piezoelectric  wafer 
active  sensors  (PWAS).  The  method  under  investigation  uses  ultrasonic  principles  where  PWAS  sensors 
are  used  to  introduce  Lamb  waves  in  the  structure  and  also  to  receive  the  reflected/transmitted  signals.  In 
this  way  the  sensors  are  used  both  as  a  transmitter  and  a  receiver.  An  array  of  PWAS  sensors  will  be  used 
to  generate  an  ultrasonic  “beam”  that  will  scan  the  area  under  investigation.  The  collected  data  will  then 
be  transferred  and  stored  into  a  PC  where  signal  processing  techniques  developed  at  USC  will  be  used  to 
identify  and  localize  the  damage.  The  key  steps  when  planning  the  POD  experiment  are  discussed  in  the 
next  sections. 


7, 1,1,1,  Establish  the  Type  and  Size  of  the  Damage  to  be  Detected 

The  specimen  under  investigation  is  a  1220x1220  mm  aluminum  plate  with  a  thickness  of  1  mm.  The 
damage  will  be  simulated  by  pin-holes  of  different  diameter.  The  holes  will  be  generated  using  a  hand- 
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held  power  drill.  We  will  eonsider  pin-holes  with  diameters  ranging  from  0.5mm,  1mm,  1.57mm  and 
2mm. 

7.LL2. 

7. 1.1.3.  Establish  the  Confidence  Interval  (Cl) 

The  probability  of  deteetion  of  a  required  damage  size  is  eorrelated  to  a  eertain  eonfidenee  interval, 
usually  90/95  %.  We  will  eonsider  in  our  experiment  a  95%  eonfidenee  interval.  In  this  way  we  want  to 
be  sure  that  we  have  a  high  deteetion  probability  and  also  a  low  false  alarm  rate. 


7. 1.1.4.  Orientation  of  the  damage 

This  experiment  will  foeus  on  two  loeations  for  the  damage.  These  loeations  are  presented  in  Figure  83. 


Figure  83  Damage  location;(a)  Phase  one,  damage  is  located  perpendicular  to  the  PWAS  array;  (b)  Phase 
two,  damage  is  located  at  an  angle 


7. 1.1. 5.  Manufacturing  the  Specimens  with  and  without  Damages 

Specimens  must  contain  a  distribution  of  damages,  as  well  as  large  areas  without  any  damages.  Because 
the  real  damaged  aircraft  specimens  are  expensive  to  duplicate,  artificial  methods  can  be  used  to  create 
damage  such  as  electrical  discharge,  machined  notches,  etc.  To  manufacture  the  damaged  specimen  we 
will  use  a  hand-held  power  drill.  The  specimen  without  the  damage  (pristine)  consists  of  an  aluminum 
plate  having  the  same  dimensions  as  the  damaged  one,  but  without  any  holes. 


7. 1.1. 6.  Hardware  and  Software  Testing 

Before  starting  the  actual  test,  it  is  important  to  completely  test  the  hardware  and  its  associated  software 
to  insure  proper  operation.  A  schematic  of  the  hardware  is  presented  in  Figure  84. 
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Figure  84  Experimental  setup  for  the  POD  experiment 


7.7.7. 7.  Documentation  of  the  Experiment 

This  step  refers  to  keeping  traek  of  all  the  methods  applied  during  the  test.  It  is  important  to  keep  traek 
and  doeument  all  the  methods  used  during  the  test.  This  step  has  a  partieular  importanee  beeause  for  an 
implementation  of  the  method  on  an  aireraft  the  approving  authority  has  to  approve  any  proeedure  that 
will  be  used  during  the  test. 


7. 1.1.8.  Environmental  Conditions 

The  POD  experiment  must  be  eondueted  in  environmental  and  physieal  eonditions  that  are  as  elose  as 
possible  to  the  real  eonditions.  In  our  initial  analysis,  we  will  use  standard  lab  eonditions. 


7. 7. 7. 9.  Calibration  Standards 

The  need  for  uniform  ealibration  standards  is  important.  By  using  the  same  ealibration  standards  the 
variables  introdueed  by  multiple  ealibration  standards  are  eliminated  and  inereases  the  level  of 
eonsisteney  with  the  results.  Not  eondueting  a  POD  involves  great  risks  that  ean  not  be  afforded.  Sueh 
risks  inelude  not  finding  a  eritieal  damage,  high  rate  of  falls  alarms,  or  using  a  wrong  teehnique  for 
damage  deteetion. 

7.1.2.  PWAS  POD  Experimental  Planning 

Deteetion:  PWAS  phased  array  using  the  EUSR  algorithm 

Pin-Hole  Damage:  Indueed  using  small  drill  bits  of  inereasing  size:  0.5  mm,  1  mm,  1.57  mm  and  2  mm 

Proeedure:  Start  with  the  smallest  drill  bit  and  inerease  until  deteetion  is  obtained  in  broadside  loeation. 
Repeat  for  offside  loeation. 
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(a) 


Figure  85  Tools  for  the  PWAS  POD  experiment:  (a)  aluminum  plate  with  a  PWAS  array  and  a  pin-hope 
damage  of  increasing  size;  (b)  set  of  small  drill  bits  in  the  size  range  0.5mm,  1mm,  1.57mm 
and  2mm 


7.2 


Results  AND  Discussions 


The  ability  to  detect  pin-hole  damage  using  PWAS  phased  arrays  and  the  EUSR  algorithm  was  tested 
during  this  procedure.  The  purpose  was  to  identify  the  minimum  pin-hole  that  can  be  detected.  The 
experiment  used  pin-hole  diameters  in  the  range  0.5  mm,  1  mm,  1.57  mm,  and  2  mm. 

The  experiment  showed  that  the  present  method  cannot  detect  pin-holes  with  0.5  mm  and  1  mm,  but  it  can 
detect  the  pin  hole  with  1.57mm  diameter.  See  Figure  86  (a).  The  tiny  dot  at  the  broadside  position  is  the 
simulation  of  the  pin  hole  in  the  2D  plot.  The  Figure  86  (c)  shows  the  result  of  the  pin-hole  with  2mm 
diameter,  which  is  much  clearer  than  the  (b). 
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d=  1.0  mm 
(/=  1.57  mm 
d  =  2.0mm 
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Figure  86  -  pin-holes:  (a)  schematic;  (b)  EUSR  GUI  mapped  image  of  the  minimum  detectable  pin-hole 
with  1 .57mm  diameter;  (c)  EUSR  GUI  mapped  image  of  pin-hole  with  2mm  diameter 
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73.  Task 4  -  Conclusions 

In  this  task  we  performed  an  initial  POD  analysis  of  the  deteetion  of  damage  using  the  PWAS  teehnology. 
A  proeedure  for  performing  the  POD  experiments  was  established.  It  was  also  shown  that  PWAS  phased- 
arrays  in  eonjunetion  with  the  EUSR  algorithm  ean  deteet  small  pin-hole  damage  as  small  as  1.57  mm  in 
a  simple  plate  strueture.  In  the  future,  this  initial  analysis  ean  be  extended  in  the  following  direetions: 

(a)  refinement  of  the  signal  proeessing  algorithm  to  bring  down  the  deteetion  limit  to  1  mm  and  even 

0.5  mm 

(b)  extension  of  the  experiment  to  eomplieated  aerospaee  struetures  with  multiple  eonstruetive  features 

(e)  eomparison  of  the  deteetion  threshold  obtained  with  the  present  method  versus  the  deteetion 
threshold  obtained  with  eonventional  ultrasonie  NDE  methods.  The  eomparison  should  inelude 
both  deteetion  threshold,  on  one  hand,  and  duration  of  test  and  manpower  requirements,  on  the 
other  hand. 
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9.  LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


1 .  PWAS  -  Piezoelectric  Wafer  Active  Sensor 

2.  EUSR  -  Embedded  Ultrasonic  Structural  Radar 

3.  SHM  -  Structural  Health  Monitoring 

4.  NDE  -  Nondestructive  Evaluation 

5 .  POD  -  Probability  of  Detection 

6.  E/M  -  Electromechanical 

7.  So  -  Symmetric  Mode 

8.  Ao  -  Antisymmetric  Mode 

9.  NDT  -  Nondestructive  Testing 

10.  POEM  -  Probability  of  false  alarm 

11.  Cl  -  Confidence  Interval 

12.  GUI  -  graphic  user  interface 


Symbols 

()i^=  imaginary  part  of  ( ) 


Orc  =  real  part  of  ( ) 


=  Fourier  Transfer  of 


() 


a  =  IJl  =  half-length  of  the  PWAS 

Ax  and  B2  =  Arbitrary  antisymmetric  motion  constants 

A2  and  Bx  =  Arbitrary  symmetric  motion  constants 

uTq  =  the  total  force  transmitted  by  the  PWAS  to  the  structure. 

c  =  chord 

c  =  Wave  Speed 

Cp  =(A  +  2ju)/ p  =  pressure  (longitudinal) 
cl  =  p/ p  =  shear  (transverse)  wave  speeds 
d  =  half  thickness 

d  =  Spacing  of  PWAS  Phased  Array 

=  Denominator/Discriminate  for  Antisymmetric  Motion 


=  Denominator/Discriminate  for  Symmetric  Motion 
E  =  Young’s  Modulus 
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=  time-harmonic  variation 
Ej^  =  signal  Energy 
fc  =  carrier  signal  of  frequency 
G{k{co),x)  =  =  wave  transfer  function 

H{x)  =  Heaviside  step  function 


1  ^ 

H(x(t))  =  —  [  =  Hilbert  transform 

TT  ^  t-T 

—00 

kHz  =  Kilohertz 

k(cD)  =  frequency-dependent  wave  number 
L  =  Arc  Length 

M  =  Number  of  PWAS  in  a  phased  array 

^4 ;  =  Numerators  for  Antisymmetric  Motion 

Na^  ;  =  Numerators  for  Symmetric  Motion 


R  =  radius  of  curvature 

Crp 

Res  @  <^  =  Residue  at  Z  =  ^ 

SQ(t)  =  a  short-duration  smoothing  window  applied  to  fc  between  0  and  tp. 
=  Signal  at  Point  P 
=  Received  Signal 
Sj  =  Transmitted  Signal 


s^{t)  =  sft)cos2/r  ft  =  the  smoothed  tone-burst  signal  generated  by  one  PWAS  sincci'  = 


sin  a  , 

- =  the 

a 


Sine  Function 


T  =  — =  tone  burst  duration  V  =  =  harmonic  electric  voltage 

/o 

u{co,x)  =  p{co)  •  G{co,x)  =  frequency  domain  response 

Ve'^"^  =  time-harmonic  voltage 
Y=  Yield  stress 

^  0.5x^  [1  -  cos(ft))]  =  Hanning  windowed  signal 
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{(j))  =  Time  of  travel  difference  for  PWAS  Signal 
5^  =  time  delay  of  m'*  PWAS  Signal 
A(x)  =  Dirac/Delta  Function 
(j)  =  phase  array  ray  angle 
=  phase  array  beam  angle 
X  =  wave  length 

\iE  =  Microstrain  (10'^  mm/mm  or  in/in) 

^  =  the  wave  number 

p  =  mass  density 

Gu  =  Ultimate  tensile  stress 

r  =  Tq  =  interfacial  shear  stress 

=  interfacial  shear  stress  between  PWAS  and  structure 

=  time-harmonic  interfacial  shear  stress 
=  Shear  Stress  in  the  y  and  x  direction 

CD  =  Angular  Frequency 
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